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I.  INTRODUCTION  AND  STATEMENT  OF  THE  PROBLEM 


The  collision  of  high  speed  projectiles  with  solid  targets  has  been  a 
subject  of  considerable  Interest  to  ollltary.  Industrial  and  academic 
communities  for  the  past  four  decaaes.  World  War  II  stimulated  analytic 
research  oriented  towards  Improving  armor  and  warhead  designs.  The  advent  of 
the  space  program  In  the  late  fifties  created  a  need  for  the  protection  of 
spacecraft  from  meteorite  Impact.  Industrial  applications  of  Impact  and 
penetration  processes  Include  mining  operations,  oil  well  stimulation,  plate, 
cable  and  pipe  cutting,  demolitions  and  many  others. 

The  subject  of  Impact  mechanics  Is  quite  broad.  It  Incorporates  such 
classical  fields  as  fluid  mechanics,  solid  mechanics,  solid  state  physics, 
mathematics  and  materials  science  Into  the  analytical  and  experimental 
treatments  of  complex,  practical  problems.  Several  general  reference  texts 
are  available  that  review  the  many  impact  type  applications  as  well  as  the 
research  techniques  developed  and  applied  with  relative  success  (Backman  [I], 
Billlngton  [2],  Goldsmith  [3],  Johnson  {4],  Klnslow  [5],  Zukas  [6],  etc.). 
Review  articles  dealing  primarily  with  spacecraft  protection  can  be  found  in 
the  Proceedings  of  the  Hypervelocity  Impact  Symposia  [7],  held  during  the 
early  sixties  (Allison  [8],  Cook  [9],  Hopkins  [10],  etc.)  and  In  the  open 
literature  (Backman  [11],  Bjork  [12],  Cook  [13],  Eichelberger  [14],  Herrmann 
[IS,  16],  Prater  [17],  Vinson  [18],  etc.).  Review  articles  concerned  with 
military  armor  applications  can  be  found  in  government  publications  (Bethe 
[19],  Perez  [20],  Wright  [21],  etc.)  as  well  as  in  the  open  literature 
(Bymslde  [22],  Christman  [23,  24],  Eichelberger  [25],  Gehrlng  [26],  Jonas 
[27],  Sagamonyan  [28],  etc.).  Both  of  the  aforementioned  categories  Involve 
the  penetration  of  high  speed  projectiles  Into  target  materials  (armor  or 
shields).  The  space  application  (conventionally  termed  hypervelocity  Impact) 
typically  Involves  geometrically  compact  projectiles  (meteorites)  with  impact 
velocities  ranging  from  1  to  80  km/s  and  targets  configured  from  thin  metallic 
plates.  Military  armor  on  the  other  hand  Is  considerably  more  massive  often 
involving  metal  target  thicknesses  on  the  order  of  hundreds  of  centimeters. 

The  projectiles,  which  these  armors  are  designed  to  stop,  usually  have  long 
cylindrical  geometries  and  Impact  velocities  significantly  less  (1-12  km/s) 
than  the  meteor  class.  In  both  categories  (hypervelocity  Impact  and  armor 
penetration)  however,  the  penetrator  velocities  are  sufficiently  great  to 
produce  Impact  stresses  much  greater  than  the  resisting  strength  of  the  target 
materials.  Consequently,  a  crater  forms  within  the  target  In  a  time  frame  of 
the  order  of  microseconds.  The  penetrator  may  disintegrate  if  the  Impact 
velocity  approaches  the  12  km/s  range;  it  will  certainly  deform  within  the 
range  of  1-12  km/s,  the  extent  of  which  will  depend  on  the  relative  magnitudes 
of  the  Impact  stresses  and  characteristic  strengths.  Hypervelocity,  within 
the  Impact  mechanics  community.  Is  usually  defined  as  a  projectile  velocity 
which  exceeds  the  longitudinal  wave  speed  of  the  target  material.  Most  metals 
have  wave  speeds  between  2  and  7  km/s.  Therefore,  depending  on  the  target, 
some  typical  Impact  velocities  may  be  less  than  the  hypervelocity  definition. 
However,  within  the  context  of  this  report,  hypervelocity  will  define  a 
velocity  sufficient  to  both  deform  the  penetrator  and  produce  a  target  cavity. 
For  most  conventional  metallic  penetrators,  an  Impact  velocity  greater  than  I 
km/s  will  be  considered  hypervelocity. 

Figure  1  depicts  three  typical  situations  Involving  the  impact  of 
representative  projectiles  upon  thin  and  thick  targets.  The  projectile 
trajectories  are  normal  to  the  target  surfaces. 
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Figure  1.  Three  Impact  Configurations,  Before  and  After  the  E/ent 

Configurations  A  and  B  represent  the  hypervelocity,  meteorite  Impact 
category.  As  noted  In  the  already  referenced  papers,  these  configurations 
have  been  extensively  studied  experimentally  (Cook  [9],  Prater  [17],  Abbott 
[29],  Holloway  [30],  Johnson  [31],  Keys  [32],  Klneke  [33],  Pond  [34],  Turpin 
[35],  Yuan  [36],  etc.),  analytically  using  engineering  models  (Davids  [37], 
Fuchs  [38],  Oplk  [39],  Rae  [40],  Rand  [41],  Rostoker  [42],  Selzew  [43], 
Sllbersteln  [44],  Scott  [45],  Taylor  [46],  Whitesides  [47],  Zald  [48],  etc.) 
and  numerically  via  finite  dlfference/flnlte  element  methods  (Bjork  [12], 
Hardage  [49],  Rlney  [50,  51,  52],  Sedgewick  [33] ,  Wilkins  [54],  etc.).  The 
problem  Is  very  complex;  wave  propagation,  unknown  dynamic  material 
properties,  phase  changes,  fracture,  time-dependent  boundary  conditions, 
nonequUlbrlxim  thermodynamics,  etc.,  all  combine  to  prevent  a  general 
analytic  solution.  Knowledge  Is  gleaned  from  the  many  analyses  of  results 
using  varying  degrees  of  approximation.  Fortunately,  most  of  the  above 
sources  of  analytic  difficulty  have  been  Investigated  either  Independently  or 
in  specific  combinations  making  the  hypervelocity  Impact  category  of  problems 
rather  well  understood.  Such  Is  not  the  case  for  the  configuration  C  of 
Figure  1. 

While  much  of  the  behavior  evident  In  the  two  other  configurations  occurs 
also  In  long  penetrator/thlck  armor  applications,  the  problem  is  considerably 
different.  Most  obvious  Is  the  shape  of  the  resulting  crater,  which  appears 
almost  hemispherical  In  B  but  becomes  more  cylindrical  as  the  length  of  the 
penetrator  Increases.  Also,  the  process  of  cavity  enlargement,  completed 
rather  quickly  In  A  and  B,  will  take  considerably  longer  when  the  length  of 
the  penetrator  increases.  Indeed,  It  Is  this  extra  time  of  Interest  that 
hinders  the  application  of  the  finite  difference  techniques  to  the  deep 
penetration  problem.  Only  a  few  numerical  studies  have  been  published  (Harlow 


[55],  Kucher  [56,  57],  etc.)  and  consequently  much  less  is  understood, 
especially  with  respect  to  individual  problem  parameter  dependence.  This  deep 
penetration  problem  has,  however,  been  of  considerable  Interest  to  the  defense 
community  for  many  years.  Protective  armors  of  tanks  and  armored  personnel 
carriers  have  in  the  past  been  designed  with  understanding  gained  from 
combined  experimental  and  engineering  model  research  approaches.  The 
experiments  (Blrkhoff  [58],  Boyle  [59],  Brooks  [60],  Cook  [13],  Elchelberger 
[61,  62],  Fugelso  [63],  Glass  [64],  Hauver  [65],  Hohler  [66],  Johnson  [67], 
Moss  [68],  Perez  [20],  Pond  [69],  Pritchard  [70],  Summers  [71],  Welhrauch  [72, 
73],  White  [74],  etc.)  identify  physical  behaviors  and  trends  that  allow  for 
subsequent  analytic  approximations.  Perhaps  the  single  most  important 
parameter  in  the  eyes  of  the  armor  designer  is  the  depth  to  which  a  projectile 
will  penetrate.  Indeed,  most  of  the  analytic  studies  to  date  have 
concentrated  on  the  penetration  rate  behavior  (the  rate  at  which  the 
penetrator/ target  Interface  moves  along  the  penetration  trajectory). 

Christman  [24]  describes  the  penetration  process  in  terms  of  four  phases: 
transient,  primary,  secondary  and  recovery.  Figure  2  shows,  schematically, 
these  four  phases  in  terms  of  pressure-time  behavior. 


Upon  initial  Impact,  shock  waves  form,  optical  flashes  are  observed  and 
the  pressures  measured  are  quite  high.  The  shocks  spread  Into  the  target 
ahead  of  the  penetrator  with  diminishing  Intensity  due  to  both  geometrical 
divergence  and  material  effects  (Hok  [75],  Rae  [40],  Torvik  [76],  etc.).  If 
the  penetrator  has  sufficient  length,  then  the  penetration  process  enters  a 
quasi-steady  phase  where  pressures,  less  than  the  transient  case,  and  along 
the  penetration  trajectory  are  almost  constant.  An  interaction  zone, 
encompassing  the  penetrator/ target  interface,  moves  with  a  constant 
characteristic  velocity  deeper  Into  the  target,  eroding  penetrator  material  as 
it  progresses.  Eventually,  the  rear  of  the  penetrator  enters  this  moving 
interaction  zone  and  upon  being  eroded,  the  process  enters  the  secondary 
phase.  The  penetrator  Is  no  longer  Important,  the  momentum  imparted  to  the 
target  by  the  first  two  phases  Is  gradually  reduced  by  irreversible  mechanisms 
until  the  static  strength  properties  of  the  target  predominate,  allowing  for 
reversible  recovery  to  occur.  This  overall  process  has  been  observed 
experimentally  (Boyle  [59],  Christman  [24],  Elchelberger  [77])  and  in  finite 
difference  calculations  (Harlow  [55],  Kucher  [57]).  The  extent  and  relative 
Importance  of  each  Individual  phase  depends  upon  both  the  penetrator  and 
target  characteristics.  For  meteorite  impact  problems,  the  primary  or 
quasi-steady  phase  never  materializes  while  for  very  long  penetrators,  this 
particular  phase  accounts  for  the  majority  of  the  penetration  process.  The 
most  serious  threat  to  ballistic  armor  involves  both  long  rods  (fired  from 
propellant  driven  guns,  .5-3  km/s)  and  very  long  metallic  jets  (2-12  km/s), 
produced  by  the  detonation  of  shaped  charge  warheads.  Under  the  presumption 
that  the  pressures  generated  greatly  exceed  the  material  strength  of  the 
target,  the  theory  for  deep  steady  penetration  was  developed  during  the  World 
War  II  era.  A  one- dimensional,  incompressible,  streamline  analogy  yielded 
estimates  for  the  penetration  rate  and  final  penetration  depth.  The 
derivation  and  subsequent  verification  (experimental  and  computational)  along 
with  a  review  of  some  of  the  assumption  relaxations  of  this  model  are 
presented  in  Appendix  A.  Many  analytic  efforts  have  considered  the 
penetration  behavior  along  the  penetration  trajectory  (Alekseevskll  [78], 
Blrkhoff  [58],  Hill  [79],  Majerus  [80],  Perez  [81],  Sagamonyan  [82],  Tate  [83, 
84],  Wright  [85]  and  others),  but  precious  few  have  considered  the  behavior 
orthogonal  to  the  penetration  path.  Wright  [21]  points  out  the  many 
limitations  associated  with  the  penetration  model  and  notes  the  need  to 
consider  the  dynamics  of  hole  growth  In  the  formulation  of  a  more  realistic 
penetration  theory.  The  growth  of  the  target  cavity,  along  directions 
orthogonal  to  the  penetration  trajectory,  deep  within  thick  targets  Is  the 
subject  of  this  report. 


The  subject  of  this  report  Involves  the  deep  penetration  of  long  rods  (or 
jets)  into  massive  targets.  Figure  3  shows  the  essential  elements  of  such  a 
problem. 
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Figure  3.  The  Deep  Penetration  of  Hypervelocity  Rods 

A  hypervelocity  penetrator,  characterized  by  Its  length  (assumed  many 
times  greater  than  Its  lateral  dimension,  density  and  velocity,  V, 

strikes  and  subsequently  penetrates  a  massive  target  characterized  by  Its 
density,  strength  properties  and  lateral  dimension.  A  cavity,  described  by 
Its  diameter,  D^,  at  specific  depths  within  the  target  forms  behind  a 

translating  Interaction  surface  moving  with  a  characteristic  velocity,  U.  It 
is  the  objective  of  this  report  to: 

a.  Develop  a  model  of  the  dynamic  growth  of  the  cavity  dimension,  D  ,  and 

Its  final  value.  ^ 

b.  Determine  appropriate  Initial  and/or  boundary  conditions  and  solve  the 
system  for  specific  problems. 

c.  Determine  the  model's  parametric  dependence,  limits  of  applicability 
and  predictive  accuracy. 

d.  Improve  the  fundamental  understanding  of  the  penetrator/ target 
Interaction. 
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II.  FORMULATION  OF  THE  PROBLEM 


The  overall  objective  of  this  report  is  to  develop  a  model  for  the  cavity 
expansion  behavior  orthogonal  to  the  penetration  path  and  subsequently  apply 
this  model  to  practical  Impact  configurations.  Parametric  dependencies  and 
model  accuracies  need  to  be  determined.  The  analysis  may  proceed  along  three 
different  approaches:  empirical  data  correlation,  approximate  analytical 
methods,  or  multidimensional  finite  difference/ finite  element  computer  codes. 
The  empirical  approach  requires  extensive  experimental  data  bases,  the 
dimensionless  groupings  of  Important  problem  parameters  are  not  always  obvious 
and  the  understanding  of  the  material  behavior  and  physical  processes  Involved 
is  not  significantly  improved  by  the  effort.'  Examples  of  this  approach  can  be 
found  in  the  review  article  by  Backman  and  Goldsmith  [11].  At  the  opposite 
extreme  is  the  numerical  approach.  The  direct  approximation  of  the  complete 
equations  of  continuum  physics,  while  allowing  for  the  treatment  of  more 
general  problems,  has  the  disadvantage  of  obscuring  Che  effects  of  individual 
problem  parameters.  Parametric  studies  require  the  repetition  of  expensive, 
time-consuming  calculations  and  often  the  resolution  is  sacrificed  in  an 
attempt  to  reduce  the  overall  cost.  The  accuracy  of  the  numerical  solutions 
depends  not  only  on  the  discretization  errors  inherent  in  the  method  but  is 
severely  dependent  on  Che  accuracy  Co  which  the  material  behavior  may  be 
described.  Zulcas,  et  al.  [6,  27]  review  the  many  contemporary  codes,  their 
respective  successes  and  limitations.  The  intermediate  approach,  approximate 
analysis,  employs  simplifying  assumptions  which  reduce  Che  general  field 
equations  to  more  analytically  tractable  forms.  The  most  important  physical 
and  material  behaviors  are  explicitly  retained  while  less  significant 
behaviors  are  neglected.  Approximations  and  simplifications  should  be 
supported  by  experimental,  theoretical  or  numerical  arguments  so  Chat  Che 
results  may  retain  reasonable  accuracy  when  compared  with  the  actual  complex 
problem. 

This  report  will  follow  the  analytical  modeling  approach  in  hopes  of 
discerning  the  most  significant  problem  parameters  and  their  influence  upon 
Che  cavity  growth  behavior. 

The  formulation  may  be  considerably  simplified  from  Che  onset  by  the 
application  of  Che  assumptions;  material  isotropy  and  axial  symmetry.  In 
addition  to  Che  obvious  analytic  advantage  gained,  these  approximations 
simplify  Che  experimental  determination  of  material  properties.  From  a 
physical  viewpoint,  both  conditions  ace  quite  appropriate  for  Che  stated 
problem.  Material  isotropy  is  conventionally  assumed  for  homogeneous, 
amorphous  and  polycrystalllne  solids.  In  the  latter  case,  isotropic 
properties  can  be  interpreted  as  averages,  taken  over  continuum  samples 
encompassing  many  randomly  oriented  grains,  of  locally  anisotropic  properties 
of  individual  grains.  Since  most  practical  armors  are  configured  with  steel 
or  aluminum  elements,  material  isotropy  will  be  assumed  throughout  Che 
remainder  of  this  report.  Subjects  involving  fiber  reinforced  materials, 
composite  penetrators  or  target  materials  that  exhibit  significant  processing 
induced  anisotropy  are,  therefore,  precluded  from  further  study.  The  second 
major  approximation,  axial  symmetry,  is  quite  appropriate  for  the  study  of 
deep  penetration.  Nearly  all  practical  penetrators  have  cylindrical  geome¬ 
tries  and  are  designed  Co  impact  a  target  end  on  (cylinder  axis  and  trajectory 
are  coincident  with  the  target  surface  normal).  For  instances  where  Che 
penetraCor  axis  Is  not  normal  to  Che  target  surface  (oblique  impact),  Che 
geometry  near  Che  entrance  surface  of  Che  target  invalidates  Che  symmetry 
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assumption,  however,  it  has  been  observed  both  experimentally  (Weihrauch  [72, 
73]}  and  computationally  (Kucher  [86])  that  for  depths  greater  than  several 
rod  diameters  below  the  oblique  entrance  surface,  the  penetration  process 
appears  axially  symmetric.  Thin  plates  and  behavior  near  oblique  free 
surfaces  will,  however,  not  be  considered.  The  lateral  boundaries  of  the 
target  will  be  assumed  to  be  either  infinite  In  dimension  or  finite  with 
cylindrical  geometry.  Actual  targets  usually  use  rectangular  plates  but  the 
lateral  dimensions  are  so  large  when  compared  with  the  final  cavity  dimension 
that  the  Infinite  approximation  Is  valid. 

Cylindrical,  Eulerlan  coordinates  are  selected  with  the  origin  located  at 
the  point  of  Initial  Impact  upon  the  target  t6p  surface.  A  cylindrical  rod 
(or  jet)  will  penetrate  the  target  along  the  axis  of  symmetry  (z  axis) 
Imparting  to  elements  of  mass  at  each  axial  plane,  an  amount  of  radial  momenta 
sufficient  to  produce  cavity  expansion.  The  penetrator  feeds  Into  a 
translating  Interaction  region  (characteristic  mushroom  shape)  where  It 
subsequently  Inverts  and  flows  parallel  to  the  expanding  cavity  surface.  As 
noted  In  the  Introduction  and  as  derived  in  Appendix  A,  the  translating 
surface  moves  through  the  target  %rlth  an  approximately  constant  rate. 
Utilization  of  pulsed  x-ray  photographic  techniques  (Weihrauch  [72,  73], 
Fugelso  [63],  and  much  unpublished  US  Army  Ballistic  Research  Laboratory 
evidence)  enables  observation  of  this  interaction  region  during  the  penetra¬ 
tion  process.  Figure  4  Is  an  example  of  such  an  observation.  The  contrast 
due  to  differences  In  material  density  allows  the  Identification  of  the  pene¬ 
trator,  cavity  and  the  interaction  regions.  While  the  fundamental  understand 
Ing  of  the  shape  and  extent  of  the  interaction  surface  Is  not  quantitatively 
developed,  the  experimental  results  qualitatively  describe  It  as  somewhat 
hemispherical  In  shape  with  dimension  greater  than  the  penetrator  diameter. 

The  experimental  results  disclose  that  not  only  does  this  surface  propagate 
according  to  the  steady  penetration  theory  but  also  the  shape  and  extent  of 
the  interaction  surface  appear  steady  as  the  penetration  process  proceeds. 

This  Is  confirmed  by  the  finite  difference  calculations  of  Harlow  [55]  and 
Kucher  [56,  57].  The  results  of  Impact  experiments  using  transparent 
materials  (glass,  PMMA,  etc.)  show  that  preceding  the  penetrator  ’’bulb"  are 
observable  waves.  If  the  penetration  rate,  U,  exceeds  the  local  sound  speed, 
c^,  of  the  target  material,  then  the  wave  resembles  a  bow  shock  associated 

with  supersonic  flow  about  a  blunt  body  of  revolution.  Figure  5  Is  an  example 
of  such  a  case  where  the  target  material  Is  PMMA  (sound  speed,  2.7  km/s)  and 
the  penetrator  Is  a  metallic  jet  moving  with  a  penetration  rate  of  4.5  km/s. 
Figure  6  Is  an  example  of  essentially  the  same  configuration  except  that  the 
target  material  Is  glass.  Here  the  characteristic  wave  speed  exceeds  the 
penetration  rate  and  the  leading  wave  is  spherical  In  shape,  appears  to  have 
originated  at  the  point  of  initial  contact  and  Is  diverging  as  It  outruns  the 
penetrator.  These  observations  are  described  In  the  experimental  papers 
(Boyle  [59],  Cook  [9],  Keys  [32],  etc.)  and  confirmed  once  again  In  the  finite 
difference  calculations  of  Hardage  [49],  Harlow  [55],  Kucher  [56,  57]  and 
others.  The  Important  point  to  note  Is  that  in  either  case  (subsonic  or 
supersonic)  the  material  ahead  of  the  penetrator  has  already  been  accelerated 
to  some  unknown  degree  by  the  time  the  penetrator/target  interface  arrives  at 
a  specific  axial  plane.  Figure  7  Is  a  sketch  that  Identifies  these  character¬ 
istics  of  the  subsonic  and  supersonic  penetration  cases  (c^  represents  the 

elastic  wave  speed).  It  will  be  assumed  that  this  behavior  exists  also  In 
optically  opaque  targets,  specifically  metals.  The  two  most  commonly  used 
target  metals  (steel  and  aluminum)  have  elastic  wave  speeds  on  the  order  of  6 
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Figure  S.  Framing  Camera  Photograph  of  a  Metallic  Jet  Penetrating 
Polymethyl  Methacrylate,  Supersonic  (U>c  )  (Courtesy  of 
George  Hauver,  USABRL) 
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km/s,  making  the  subsonic  case  likely  for  most  practical  applications.  The 
nature  and  behavior  of  these  Impact  generated  waves  in  solids  has  been  the 
subject  of  considerable  research  and  falls  within  the  general  category  of  wave 
propagation  (Backman  [87],  Crlstescu  [88],  Davids  [37],  von  Karman  [89],  Koehler 
[90],  Kolsky  [91],  Llpkln  [92],  Mok  [75],  Wilkins  [54],  etc.}.  The  subject 
will  be  addressed  later  in  this  text  when  boundary  conditions  are  considered. 

Up  to  this  point,  little  has  been  said  about  the  nature  of  the  material 
flow  within  the  target  that  ultimately  determines  the  crater  surface. 

Experiments  utilizing  wax  and  plasticine  (Frazier  [93],  Holloway  [30],  Johnson 
[67])  and  those  using  metals  (Dong  [94],  Glass  [64],  Pond  [69],  Welhrauch  [72, 
73])  have  approached  this  topic  by  studying  postmortem  deformation  fields. 

The  embedding  of  material  tracers  within  the  target  or  the  layering  of 
dissimilar  plates  allows  the  measurement  of  displacements  and  subsequent 
calculations  of  finite  strain  tensors  and  invariants  to  be  performed.  The 
observations  from  these  experiments  may  be  siimmarlzed  as  follows: 

a.  Target  mass  remains  essentially  constant  throughout  the  penetration 
process.  Some  mass  is  ejected  during  the  initial  transient  phase  as  a  shear 
lip  forms  but  as  the  depth  of  the  cavity  Increases,  this  fraction  of  the  total 
mass  becomes  negligible. 

b.  The  target  density,  before  and  after  penetration,  varies  only  slightly 
(on  the  order  of  one  percent  or  less). 

c.  The  target  material  between  the  cavity  surface  and  a  radial  surface 
slightly  larger  undergoes  finite  radial  and  axial  displacements  (shear 
deformation) . 

d.  Beyond  this  unspecified  radius,  the  deformation  appears  to  be  pure 
radial  expansion.  The  degree  of  deformation  decreases  with  Increasing  radial 
coordinate. 

e.  For  certain  polycrystalline  targets  it  is  sometimes  possible  to 
discern  a  maximum  radius  of  permanent  grain  reorganization  (limit  of  plastic 
deformation) . 

f.  The  deformation  field  is  symmetric  with  respect  to  the  penetration 
axis. 

The  above  observations  Involve  only  the  final  deformation  field.  However, 
Welhrauch  [72],  using  a  target  composed  of  copper  plates,  was  able  to 
sequentially  observe  the  penetration  of  a  copper  rod  during  the  event.  Figure 
8A  is  a  sketch  of  these  observations.  The  surfaces  of  the  individual  plates 
appear  to  deform  in  a  parallel  fashion.  Once  again,  the  finite  difference 
calculations  of  Kucher  [56],  Figure  8B,  agree  qualitatively  well  with  these 
experiments. 


Deformation  Fields  During  the  Penetration  Event 


A.  Breakdown  into  Regions  of  Different  Deforaatlons 

With  the  support  of  both  experimental  and  numerical  evidence,  the 
following  situation  will  be  assumed:  the  deformation  of  a  target  being 
penetrated  by  a  long,  axially  aligned  rod,  separates  into  four  distinct  but 
interacting  spatial  regions.  Figure  9  Identifies  these  four  regions,  ordered 
according  to  their  relative  radial  dimensions.  The  regions  are  separated  by 
boundaries  that  translate  in  time  along  axial  planes  (z  -  constant)  as  the 
cavity  expands.  The  wave  that  precedes  the  penetrator  "bulb'*  is  described  as 
an  elastlc/plastlc  wave  which  propagates  with  a  characteristic  wave  speed.  It 
separates  quiescent  target  material  from  that  which  has  been  affected  in 
varying  degrees  by  the  passing  of  this  wave.  The  stress  distribution  within 
this  affected  zone,  although  not  quantitatively  defined,  is  expected  to  vary 
from  negligible  magnitudes  at  the  larger  radii  to  the  stagnation  pressure  at 
the  symmetry  axis.  At  some  radius,  b,  the  stresses  should  be  sufficient  to 
meet  or  exceed  standard  plastic  yield  criteria.  Region  IV  then  represents  the 
part  of  the  target,  upon  a  specific  axial  plane,  that  can  be  appropriately 
modeled  with  elasticity  theory.  Region  III,  bounded  by  the  radii  a'  and  b, 
represents  target  material  that  undergoes  plane  radial  deformations  under 
stresses  greater  than  the  elastic  range.  Region  II  Involves  plastic 
deformation  in  both  radial  and  axial  directions.  It  is  bounded  by  the  radius, 
a',  and  the  cavity  wall,  a.  Between  the  cavity  and  the  axis  of  symmetry. 
Region  I,  the  penetrator  enters  the  translating  Interaction  region,  inverts 
and  flows  along  the  target  cavity.  In  each  individual  region,  the  nature  of 
the  deformations  %flll  determine  the  appropriate  formulations.  The  four 
regions  will  be  ultimately  coupled  by  the  requirement  of  continuity  at  each  of 
the  boundaries.  For  any  speclHed  axial  coordinate,  the  cavity  expansion 
process  begins  when  the  penetrator/ target  "bulb"  passes  this  axial  plane  and 
ends  when  the  stresses  in  the  target  decay  to  the  point  where  material 
strength  halts  further  deformation. 


B.  The  Elastic  Region  IV 


From  the  principle  of  linear  momentum,  the  equations  of  motion  in 
differential  tensor  form  may  be  written: 
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where  U.  represents  the  components  of  the  velocity  vector,  v..  is  the  stress 
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tensor,  b^  is  the  body  force  vector,  p  is  the  density  and  the  operator  is 

the  material  derivative  operator.  Cylindrical  Eulerlan  coordinates  have 
already  been  selected  and  with  the  assumption  of  plane  axlsymmetrlc 
deformation,  the  radial  equation  of  motion  is: 
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where  a  and  <r..  are  the  radial  and  circumferential  stresses,  respectively, 
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which  are  also  principal  in  this  configuration.  If  inertia  forces  are  assumed 
negligible  throughout  this  region  (they  certainly  are  smaller  than  those  in 
Regions  II  and  III),  then  in  absence  of  applied  body  forces,  equation  (2) 
simplifies  to  the  plane  strain  equation  of  equilibrium: 


Figure  9.  Sketch  of  the  Separation  of  the  Target  Into  Regions 
of  Different  Deformations  (Subsonic  Flow) 
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which  is  appropriate  for  the  elastic  loading  of  a  long,  thick  tube  by  uniform 
internal  pressure.  The  solution  to  this  problem  is  well  known  (Timoshenko 
[95]).  The  stress  distribution  is  of  the  form: 
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where  v  is  Poisson's  ratio,  and  are  constants  determined  by  application 
of  the  boundary  conditions  and  the  axial  stress,  «  ,  is  not  zero  in  general. 

2Z 

If  the  target  has  relatively  large  lateral  dimension  such  that  the 
semi-infinite  approximation  is  valid,  then  it  may  be  assumed  that  the  radial 
stress  vanishes  for  large  radial  coordinates.  The  constant,  K^,  must 

therefore  be  zero  and  the  stress  distribution  for  the  semi-infinite  target  is 
given  by: 
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If  the  lateral  boundary  has  a  finite  dimension,  c,  the  condition  that  the 
radial  stress  vanishes  at  this  boundary  yields  a  slightly  different  form: 

Orr  -  Ko  <1  - 
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The  remaining  undetermined  constant,  for  either  case,  requires  that 
another  boundary  condition  be  specified.  This  will  be  accomplished  once  the 
formulation  in  Region  III  is  completed  and  stress  continuity  at  the 
elastic/plastic  boundary,  b,  is  imposed. 

C.  The  Plane  Plastic  Region  III 

The  elastic  Region  IV  utilized  the  assumptions  of  axisymmetrlc  plane 
deformation,  negligible  inertia  effects  and  small  displacements  in  order  to 
determine  the  internal  stress  distributions.  The  last  two  approximations  are 
perhaps  only  accurate  for  large  radial  positions,  sufficiently  far  from  the 
penetrator/ target  Interaction  "bulb.”  Region  III  encompasses  material  of 
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significantly  smaller  radii,  thereby  precluding  the  use  of  the  same  set  of 
assumptions.  Experimental  evidence  supports  the  contention  of  a  region  of 
axisymmetrlc  planar  flow,  but  also  categorizes  the  deformation  as  finite  and 
irreversible.  The  experiments  additionally  observe  that  the  flow  process  does 
not  include  substantial  permanent  density  change.  The  numerical  studies  show 
that  compressibility  effects,  while  not  completely  negligible,  are  usually 
spacially  localized  around  shock  surfaces  which  eventually  diverge  with  time 
and  distance.  It  Is  expected  that  strains  as  well  as  strain  rates  will  be 
large  within  this  region  but  It  will  be  assumed  that  throughout  the  process, 
incompressibility  will  apply.  The  analysis  will  be  limited  to  plastic,  radial 
deformations. 


Let  (r,  z,  $,  t)  be  the  Instantaneous  velocity  components  of 

plastically  flowing  target  material  at  any  given  Eulerlan  coordinate  point  (r, 
6,  z) .  The  axial  symmetry  of  the  problem  requires  that  the  tangential 
velocity  component,  U^,  vanish  and  the  presumption  of  plane  flow  eliminates 

the  axial  component,  U  .  Then  with  respect  to  a  specific  axial  plane,  only 

z 

the  radial  velocity,  U  ,  remains.  The  continuity  equation  for  an 
incompressible  flow  Is: 


(7) 


where  U  varies  with  radial  position,  r,  and  time.  The  solution  to  this 
equation  Is: 
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where  0(t}  represents  an  arbitrary  but  differentiable  function  of  time.  The 
tenser  equations  of  motion,  in  terms  of  the  stress  tensor,  ,  the  scalar 

density,  p  ,  the  coordinates,  x^ ,  velocity  vector,  ,  and  body  force  vector, 

b^,  may  be  written: 
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In  the  absence  of  body  forces,  the  radial  equation  of  motion  for  plane 
axisymmetrlc  flow  simplifies  with  <7rz>  '^z  ~ 
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As  in  the  plane  elastic  Region  IV,  the  stresses  are  also  principal  In  this 
configuration.  Substitution  of  equation  (8)  into  equation  (10)  yields: 
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where  the  dots  denote  differentiation  with  respect  to  time  and  the  stress  is 
described  by  the  principle  stress  vector.  Equation  (11)  is  identical  to  the 
governing  equation  of  motion  derived  by  Hill,  et  al.  [79]  and  later  by  Appelby 
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[96]  and  Al-Hassanl  [97]  for  a  thick  walled.  Incompressible  axially 
constrained  tube  under  uniform  Internal  pressure.  Appelby's  formulation 
assumed  plane  strain  deformation,  stress  free  outer  boundaries,  Tresca  yield 
behavior  and  a  vlsco-plastlc  flow  rule.  Expansion  histories  of  the  Inner 
surface  were  prescribed  and  the  internal  pressure,  necessary  to  produce  such, 
was  calculated.  A  comparison  was  made  of  the  Individual  contributions  of 
Inertia,  viscosity  and  perfect  plasticity  to  this  time  dependent  pressure. 

The  Inverse  problem  of  specifying  the  pressure  history  and  then  solving  for 
the  velocity  and  deformation  fields  was  not  attempted.  An  earlier  Russian 
effort,  Agagablan  [98],  considered  the  same  geometrical  configuration  and  very 
similar  relations  for  plastic  yield  and  flow  but  allowed  only  small 
displacements.  The  radial  equation  of  motion  derived  by  this  author  was: 

(12, 
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which  differs  from  equation  (11)  in  the  acceleration  term,  lacking  the 
contribution  of  the  convective  term.  With  Identical  boundary  conditions  as 
Appelby,  Agagablan  obtained  an  exact  solution  for  the  stress  components 
(spatial  and  temporal).  Agagablan  also  considered  the  case  of  a  very  thick 
tube  such  that  the  entire  tube  was  not  undergoing  plastic  deformation.  The 
boundary  separating  elastic  and  plastic  deformations  was  assumed  continuous 
with  respect  to  displacements  and  radial  stress.  The  pressure  was  applied  In 
a  Heaviside  step  fashion  with  arbitrary  magnitude.  Stress  distributions  as 
well  as  the  propagation  behavior  of  the  elastlc/plastlc  Interface  were 
analytically  determined.  However,  the  chief  limitations  of  this  analysis, 
when  considering  the  cavity  expansion  problem,  are  the  small  displacement 
approximation  and  the  specification  of  unknown  pressure  history.  Since  this 
report  Is  concerned  with  large  radial  displacements,  the  formulation  will 
follow  an  approach  parallel  td.th  the  work  of  Hill,  Appelby  and  Al-Hassani . 

Returning  to  equation  (11),  the  term  furthest  to  the  right  hand  side 

Involves  the  difference  of  the  two  principle  stresses,  v  and  a..  The 

r  tr 

Tresca-Hohr  yield  criterion  (condition  of  constant  maximum  shear  stress) 
presumes  that  the  greatest  difference  between  the  principle  stress  components 
equals  twice  the  Inherent  material  shear  strength,  r,  (or  the  uniaxial  yield 
strength  In  simple  tension,  If  the  body  Is  deforming  plastically.  Within 

the  plane  deformation  regions  of  this  problem,  this  criterion  may  be  stated 
as: 
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where  all  of  the  Inequalities  hold  for  the  elastic  state  but  only  one  or  two 
equalities  hold  while  the  material  is  in  the  yield  state.  The  selection  of 
the  Tresca-Mohr  criterion  over  the  commonly  used  Mlses-Hencky  criterion  Is 
typically  based  upon  the  ease  to  which  either  criterion  may  be  applied  to  a 
specific  problem.  Both  are  usually  regarded  as  equally  valid  and  the 


differences  observed  are,  for  Che  most  part,  Insignificant.  Over  the  past 
years,  Che  Tresca-Mohr  hypothesis  has  been  found  to  agree  reasonably  well  with 
experimental  evidence  and  to  satisfactorily  characterize  the  yield  state.  Due 
to  Its  conceptual  simplicity  and  ease  of  application  to  this  problem,  Che 
Tresca-Mohr  condition  will  be  assumed  throughout  the  remainder  of  this  report. 

It  is  necessary  to  determine  which  equality  (13-15)  Is  applicable  for 
target  material  In  the  yield  state.  The  maximum  shear  stress  is  half  of  Che 
difference  of  the  greatest  and  least  principal  stresses;  the  Intermediate 
principal  stress  does  not  Influence  the  yield  state.  An  ordering  of  the  three 
principal  stresses,  therefore,  determines  the  appropriate  equality  (13-15). 
Reference  to  Che  equilibrium  stress  distributions  of  thick-walled  Cubes 
subject  to  internal  pressures  sufficient  to  produce  plastic  deformation  (Allen 
and  Sopwlth  [99],  Bethe  [19],  Cook  [100],  Hill,  Lee  and  Tupper  [101],  Kachanov 
[102],  Nadal  [103],  etc.)  provides,  perhaps.  Insight  for  Che  dynamic  problem. 
In  all  of  these  references,  Che  radial  and  circumferential  stress  distribu¬ 
tions  are  essentially  Che  same  but  the  axial  stresses  differ  according  to  the 
assumptions  chosen  to  allow  solution.  Kachanov  notes  that  the  axial  stress  is 
half  Che  sum  of  Che  radial  and  circumferential  principal  stresses  for  both  the 
elastic  state  and  Che  pure  plastic  state  of  a  thin-walled  tube.  Nadal  and 
subsequently  Cook  proposed  Chat  Indeed  this  relationship  should  be  a  very  good 
approximation  for  Che  plastic  state  of  a  thick-walled  cylinder.  Hill,  et  al., 
using  a  finite  difference  technique  to  solve  Che  boundary  value  problem  of 
pressure  applied  on  the  inner  surface  and  a  traction  free  exterior  boundary, 
found  Chat  as  plastic  flow  develops  (in  terms  of  increasing  deformation),  Che 
axial  stress  approaches  Che  mean  of  the  two  ocher  stresses,  especially  for 
elements  closest  Co  the  Inside  Cube  surface.  Allen  and  Sopwith,  after 
reviewing  the  solutions  available  at  their  time,  compared  the  independent 
predictions  with  particular  reference  to  the  exterior  boundary  expansion  and 
concluded  that  Cook's  approximation  of  the  axial  stress  relationship  leads  to 
only  Insignificant  errors.  Additionally,  Bethe  pointed  out  Chat  plastic  flow 
in  Che  outer  regions  has  already  been  assumed  to  be  axially  independent,  the 

greatest  stress  difference  Is,  therefore,  l<r  'v.l.  This  Is  consistent  with 

r  cf 

Cook's  assumption  and  has  been  qualitatively  confirmed  for  Che  dynamic  case  by 
Che  finite  difference  calculations  of  Wagner  [104].  The  equality  (13)  will 
Chen  represent  Che  condition  of  plastic  flow: 

“  °y 

It  would  be  quite  convenient  at  this  point  that  Che  yield  strength,  ,  be 
a  material  property  independent  of  the  strain  history,  temperature,  ^ 

hydrostatic  pressure,  etc.  Unfortunately,  this  Is  not  Che  case  for  most 
situations  Involving  Impact  dynamics.  Appendix  B  gives  a  brief  review  of  some 
experiments  and  modeling  efforts  concerned  with  the  variation  of  this  yield 
strength  during  Impact  Initiated  deformation.  One  of  Che  most  obvious 
dependencies  of  Che  yield  strength  is  upon  the  rate  of  strain  (as  several 
models  reflect).  Perhaps  one  of  Che  more  widely  utilized  models,  especially 
within  Che  Russian  literature  (Agagablan  [98],  Godunov  [105,  106],  etc.)  is 
Che  vlscoplastlc  flow  law  originally  proposed  by  Hohenemser  and  Prager  [107]. 
It  can  be  written  in  the  following  forms  (ex.  Prager  [108]): 


(17) 
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where  e^^  represents  the  components  of  the  plastic  strain  rate  tensor,  (r^j  1: 

the  stress  tensor,  n  Is  the  constant  coefficient  of  viscosity  and  F  Is  the 
yield  function  defined  as: 


F  -  |oi-a2|  - 


and 


<F> 

<F> 


0  if  F  <  0 
F  if  F  >  0 


(18) 


where  and  maximum  and  minimum  principal  stress  components 

respectively.  Following  the  convention  that  tensile  stresses  are  positive  and 
recognizing  from  the  equilibrium  solutions  that  ,  the  yield  function  (18) 

V  r 

and  equation  (16)  yield: 


F  -  09-Or-Oy 


(19) 


with  this  specific  yield  function,  equation  (17)  reduces  to  two  nonvanishing 
components  of  plastic  strain  rate: 

ueJr  "  Or‘<J9+Oy  (20^ 

and 
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As  Is  typically  assumed  in  Incremental  plasticity  (ex.  Fung  [109]),  the 
rate  of  deformation  tensor  may  be  decomposed  Into  elastic  and  plastic 
components : 
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and  for  most  problems  Involving  large  plastic  deformations,  which  this 
report  subject  Is  certainly  one,  the  elastic  contribution  is  assumed 
negligible,  therefore. 


*ij 


The  nonzero  components  of  the  rate  of  deformation  tensor  for  plane 
axlsymmetrlc  flow  can  be  written  as: 


(23) 


Differentiating  and  substituting  equation  (8)  Into  equations  (24)  and  (25)  and 
combining  the  results  with  equations  (23),  (20)  and  (21)  yields  the  following 
relation: 

U  ^2  °e-<JT-Oy  •  (26) 

Substitution  of  this  equation  Into  equation  (11)  then  produces  an  expression 
relating  the  radial  stress  gradient  and  the  function  0  (t): 


which  may  be  directly  Integrated  with  respect  to  r,  between  the  limits  of  the 
radial  coordinate,  r'  (a*  £  r'  £  b)  and  the  elastic/plastic  boundary,  b. 


r(r')-Or(b)  -  j0(»+Oy|  \  ^ 

The  circumferential  stress,  <r  ,  from  equation  (26)  Is  similarly  written: 

)  -  Oy+Or(*>)  +  jp0+Oy|  ^  +  -I  p  "  ^2)  ^  ^ 


O0(r')  -  Oy+Cr^h)  +  jp0+OyJ  —  +  2  ^  ^ 

The  Instantaneous  stress  distribution  Is  thereby  determined  for  time,  t',  and 

any  position  within  Region  III  If  the  instantaneous  position  of  b(t')  and  the 

values  0(t'),  Cl(t'),  d(t')  are  known.  At  this  point,  0(t)  remains 
undetermined. 

D.  The  Three  Dimensional  Plastic  Region  II 

In  the  Introductory  portion  of  this  chapter,  observations  of  final  target 
deformations  from  experimental  and  numerical  investigations  were  summarized. 
This  region  of  the  target,  nearest  to  the  cavity  surface,  exhibits  axial  and 
radial  deformation.  The  degree  of  deformation  Is  great  and  It  Is  also 
expected  that  the  deformation  rates  are  large.  Since  the  Irreversible 
deformations  are  so  much  greater  than  elastic  magnitude,  plastic  behavior 
should  dominate.  Glass  [64]  and  Pond  [69]  experimentally  determined  the  first 
strain  Invariant  within  the  region  adjacent  to  cavities  In  steel  targets. 

Even  at  the  cavity  surface,  where  strains  are  maximum,  the  final  compression 
was  less  than  one  percent.  Similar  observations  were  made  by  Holloway  [30] 
who  considered  a  much  more  compressible  material,  wax.  Therefore,  In  addition 
to  the  approximations  of  axial  symmetry  and  material  Isotropy,  the  plastically 
deforming  target  material  In  Region  II  will  be  assumed  Incompressible.  Since 
the  deformations  are  not  restricted  to  axial  planes,  shear  stress  and  shear 
flow  are  likely  to  some  unknown  degree  in  the  r-z  plane. 

The  differential  equations  of  motion  for  axlsymmetrlc  flow  in  Eulerlan 
coordinates  may  be  written: 
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where  u  ,  u.,  and  u  are  the  Instantaneous  components  of  velocity,  p  is  the 
Xu  z 


target  density  and  are  the  nonzero  components  of  the  instantaneous  stress 


tensor.  The  incompressible  continuity  equation  is: 

-  0  . 


3u-  u.  3u- 

-T^  +  —  + 

3r  r 


3z 


(33) 


The  system  of  nonlinear  partial  differential  equations  (30  -  33)  is 
indeterminate.  Either  additional  Independent  equations  must  be  supplied  or 
assumptions  that  reduce  the  number  of  dependent  variables  must  be  Imposed  in 
order  to  obtain  solutions.  Most  contemporary  finite  difference  approaches 
utilize  the  former  method.  The  incompressibility  condition  is  relaxed,  an 
energy  conservation  equation  is  formulated,  constitutive  equations  are  assumed 
and  the  resulting  more  complex  system  of  nonlinear  partial  differential 
equations  are  approximated  with  a  set  of  algebraic  difference  equations  which 
are  then  solved  on  the  computer  in  a  variety  of  ways.  Alder  [110]  describes 
some  of  these  and  presents  examples  of  their  relative  accuracies. 
Unfortunately,  this  general  approach  is  not  without  difficulties.  In  addition 
to  the  errors  introduced  by  the  differencing  approximation,  the  formulation  of 
the  energy  equation,  the  choice  of  the  most  appropriate  constitutive  model  and 
the  treatment  of  boundary  conditions  and  shock  singularities  all  involve 
varying  degrees  of  approximation.  (k)st  is  often  a  limitation,  especially  when 
high  resolution  and  minimal  errors  are  desired.  The  alternate  approach  to  the 
solution  of  the  system  of  field  equations  is  to  reduce  the  original  set 
(30-33)  to  simpler  analytic  form  by  suitable  approximations.  These 
simplifications  should  be  supported  by  experimental  evidence  or  theoretical 
argument . 


The  experiments  of  Glass  [64],  Moss  [111]  and  Pond  [34,  69]  utilized  some 
metallographlc  techniques  to  study  the  deformation  fields  of  the  target  near 
hypervelocity  Impact  cavities.  While  Pond  and  Glass  embedded  markers  In  the 
undeformed  target  to  trace  point  displacements.  Moss  recognized  chat  the 
microstrucCure  typical  of  quenched  and  tempered,  hardenable,  low-carbon  steel 
plate  could  serve  as  a  reference  grid  from  which  shear  strains  can  be 
calculated.  Bands  of  plane  chemical  inhomogeneicies  Chat  had  been  spread 
through  Che  plate  as  it  was  rolled  from  an  ingot  are  initially  parallel  with 
Che  Cop  and  bottom  surfaces  of  Che  target  plate.  These  reference  bands  are 
deformed  as  the  penetraCor  perforates  a  given  axial  plane.  Measurements  of 
Che  slopes  of  the  reference  bands  then  determine  the  shear  strain 
distribution.  Moss  investigated  Che  nature  of  adiabatic  shear  bands  by  shear 
plugging  plates  with  explosively  driven  punches.  He  estimates  the  shear 
strain  rates  and  temperatures  within  these  extremely  small  regions  (on  Che 
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order  of  10  cm)  which  run  orthogonal  to  the  metallographlc  reference  bands. 
Immediately  adjacent  to  these  shear  bands  is  a  region  which  resembles  Che 
deformation  fields  observed  by  Welhrauch  [72]  and  Johnson  [67]  and  calculated 
by  Kucher  [56] .  The  process  of  shear  plugging  is  considerably  different  from 
long  rod  penetration.  Consequently,  an  experiment  was  proposed  and  performed 


utilizing  the  metallographic  technique  of  Moss,  but  applied  to  the  steady 
deep  penetration  problem.  Figure  10  shows  the  experimental  configuration 
designed  to  exhibit  the  deformation  field  in  Region  II. 


Figure  10.  Test  Configuration  for  the  Observation  of  the  Final 

Deformation  Field 

The  detonation  of  a  shaped  charge  device  produces  an  axially  aligned, 
continuous,  metallic  jet  with  a  velocity  range  typically  between  2-8  km/s. 

The  top  target  plate  was  located  close  enough  to  the  warhead  so  that  the  Jet 
penetrated  the  target  in  a  continuous  state  (at  great  distances  the  jet 
particulates  due  to  the  Inherent  axial  velocity  gradient).  The  warhead  was 
selected  so  that  the  target  would  be  completely  penetrated.  According  to 
Christman  [23],  upon  Initial  Impact  with  the  top  surface,  the  Jet  begins  the 
transient  penetration  phase  which  quickly  disappears  and  then  enters  the 
primary  or  quasi-steady  penetration  mode  as  described  In  Appendix  A.  After 
the  event  was  completed,  the  perforated  plates  were  recovered  and  Inspected. 
Plate  number  two  was  selected  for  further  study  since  It  Involved  only  the 
quasi-steady  phase  and  the  jet  would  still  be  completely  continuous  at  this 
distance.  This  plate  was  then  sawed  in  two  along  the  axis  of  the  hole,  a 
section  (see  Figure  10)  was  removed,  polished  and  macroetched  with 
Oberhoeffer' s  etch.  Figure  11  Includes  a  sketch  and  photograph  of  the  exposed 
banded  structure  after  the  penetration  event,  as  observed  under  the 
microscope.  The  fibrous-llke  banded  structure  of  alternating  layers  of 
ferrite  and  pearllte  originally  aligned  along  axial  planes,  deforms  into 
directions  parallel  to  the  penetration  trajectory  near  the  cavity  surface. 

The  surfaces  appear  practically  parallel  to  each  other,  especially  as  the 
distance  from  the  cavity  Increases.  The  shear  strain,  y  ,  is  defined  as: 


(34) 


1  /3“z  ^  3“z\ 
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where  and  represent  the  final  displacements.  If  stretching  In  the  axial 

direction  Is  neglected  (pure  shear  deformation  In  the  r-z  plane)  then  the 
axial  gradient  vanishes  and  equation  (34)  reduces  to  a  relationship  between 
the  shear  strain  and  the  measured  deflection  angle, 

1  3^2  1  (35) 
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Several  reference  bands  were  selected  for  different  axial  planes  through 
the  thickness  of  plate  number  two.  Deflection  angles  were  measured  along  each 
of  these  bands  for  radial  positions  close  to  the  cavity  surface.  Figure  12 
gives  the  calculated  shear  strain  distribution  with  respect  to  the  radial 
coordinate.  All  of  the  reference  bands  yield  distributions  that  fall  within 
the  shaded  region,  decay  to  negligible  magnitudes  at  about  the  same  radial 
coordinate,  and  exhibit  no  consistent  variation  with  respect  to  axial 
position. 

The  observation  of  these  reference  band  surfaces  as  well  as  the 
experimental  evidence  of  Welhrauch  [72]  and  the  numerical  calculations  of 
Kucher  [56]  all  suggest  the  assumption  of  parallel  flow.  The  axial  gradients 
of  displacements  and  velocities,  therefore,  vanish  In  Region  II.  The 
approximation  does  not  require  that  either  the  shear  stress,  „  ,  or  shear 
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strain,  v  equal  zero,  al .hough  they  vanish  by  definition  In  Regions  III  and 
&  z 

IV.  Equations  (30-33)  simplify  to: 
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Equation  (38)  Is  Identical  to  the  Incompressible  continuity  equation 
applicable  to  the  plane  deformation  Region  III.  The  radial  velocity  is, 
therefore: 


but  the  parallel  flow  restriction  eliminates  the  axial  dependence  and  the 
radial  velocity  then  has  identical  form  as  In  Region  III: 


(40) 
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where  d(t)  is  an  arbitrary,  differentiable  function  of  time.  Equation  (36) 


includes  the  axial  gradient  of  the  shear  stress,  „ 


This  shear  stress  is 


given  for  the  case  of  a  Newtonian  flow  (which  is  perhaps  appropriate  in  this 
region  where  deformations  and  rates  are  greatest): 


where  m  is  the  viscosity  coefficient  as  described  in  the  formulation  of  the 
vlscoplastlc  flow  of  Region  III.  The  assumption  of  parallel  flow  once  again 
eliminates  the  axial  gradient  of  the  radial  velocity  and  upon  differentiation 
with  respect  to  z,  equation  (41)  becomes: 

The  axial  velocity  field  Is  assumed  to  be  piece-wise  continuous,  thereby 
allowing  the  interchange  of  the  partial  derivatives  and  therefore: 


Combining  equations  (36),  (37),  (40),  (41)  and  (43)  yields: 
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where  because  the  shear  stress,  «  ,  is  not  necessarily  zero,  the  stresses 


»rr* 


and  <r  are  not  principal.  The  axial  velocity  distribution  (U.(r,t)) 
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Is  undefined  at  this  point,  but  does  not  appear  explicitly  in  equation  (44). 
With  the  exception  that  the  stresses  are  not  exactly  principal,  equation  (44) 
has  the  identical  form  of  the  radial  equation  of  motion  developed  in  Region 
III.  The  stresses  in  Region  III  are  principally  due  to  the  assumption  of 
plane  deformation.  The  Tresca  yield  condition  and  the  Bingham-like  flow  rule 
of  Hohenemser  and  Prager  [107]  were  selected  because  of  their  analytic 
simplicity  and  because  they  accurately  model  behavior  of  metal  deformations 
typical  of  Impact  and  penetration  problems  (see  Appendix  B).  In  order  to 
continue  the  formulation  of  Region  II  along  the  same  approach  as  in  Region 
III,  the  relative  magnitudes  of  the  stress  components  must  be  discerned. 

Equation  (41)  suggests  that  the  plastically  flowing  target  material 
behaves  much  like  a  viscous  fluid.  The  solid  nature  of  the  target  has  been 
implicitly  assumed  negligible  in  comparison  with  the  viscous  stresses  that 
develop  because  of  the  very  large  deformation  rates  within  this  region. 
Suppose  that  the  stress  tensor  may  be  described  in  more  general  form  by: 


(46) 
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where  /i  and  represent  the  conventional  fluid  viscosity  coefficients, 


represents  the  rate  independent  components  of  the  solid  stress  tensor  and  x. 


are  the  cylindrical  coordinates.  Incompressibility  (0Uj^/8Xj^  ~  0),  axial 

symmetryf^,  U.  *  o)  and  parallel  flow  (8U  /8  >  8U  /8  ■  0)  simplify  equation 

(46)  to: 
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Utilizing  equation  (40),  the  nonzero  components  can  be  written: 


rr 


-  2u 
rr  r 


'rz 


O  +  u 
rz 


3u2 

3r 


o  «  0 

<^00  -  <^90  +  2u  ^ 


zz 


(47) 


(48) 


and  further  assuming  that  the  rate  dependent  contributions  predominate  over 
the  static  terms  during  the  time  when  rates  of  deformation  are  great,  the 
components  are  approximately  given  as: 


(0-z  planes) , 
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Depending  on  the  magnitude  of  the  radial  gradient  of  the  axial  velocity, 
the  maximum  shear  stress  appears  upon  either  axial  (r-9)  planes  or  planes 
Including  the  symmetry  axis  (r-z).  If  It  Is  assumed  that  the  shear  strain  at 
any  coordinate  develops  continuously  over  the  time  Interval  dt  during  which 
Incremental  deformations  occur,  the  shear  strain  rate  Is  then  related  to  the 
shear  strain  by: 
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or  In  other  words,  the  velocity  and  deformation  gradients  possess  Identical 
spatial  functional  forms.  Figure  12  exhibits  an  observed  shear  strain 
distribution.  The  shear  strain  rate  may  be  approximated  with: 


-  K  exp  (r-a)| 


(54) 


where  K,  the  maximum  possible  shear  rate,  can  be  estimated  by  the  presumption 
of  slip  at  the  target  cavity  when  the  shear  strength  of  the  target  Is 
exceeded. 
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The  term  it/r  represents  u^/r.  Within  Region  II,  the  radii  are  small 

(relative  to  the  target  dimension)  and  the  radial  velocities  are  maximum. 

2 

Therefore,  It  will  be  assumed  that  throughout  Region  II  that  0/r  >  and 

Assuming  this  to  be  the  case,  then  the  greatest  shear  stress 
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Is  likely  upon  axial  (r-9)  planes.  If,  however,  the  radial  gradient  of  axial 

2 

velocity  Is  of  the  same  order  as  0/r  then,  shear  deformation  Is  equally 
likely  upon  both  (r-z)  and  (r-tf)  planes. 
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The  Tresca  yield  criterion  was  applied  within  the  formulation  of  Region 
III.  The  criterion  considers  the  differences  of  the  principal  stresses,  and 
as  already  noted,  the  hoop  and  radial  stresses  In  Region  II  are  not 
necessarily  principal.  They  are  principal  If  the  radial  gradient  of  axial 
velocity  is  completely  negligible.  However,  the  results  depicted  in  Figure  12 
and  the  argument  leading  to  equation  (54)  suggest  that  this  Is  not  the  case. 

.  2 

If  the  gradient  term  is  of  the  same  order  as  0/r  everywhere  within  this 
region  (actually  this  Is  perhaps  a  worst  case  situation)  then  the  principal 
stress  state  will  be  different  In  general  from  the  actual  stress  state.  As  a 
measure  of  how  different,  let 
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The  approximate  viscous  stress  state,  equations  (49)  with  equation  (56)  are 
related  to  the  principal  stress  state,  S^,  by  the  follovirlng  relationship: 
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(57) 


The  solution  to  the  resulting  cubic,  algebraic  equation  yields  three  real, 
distinct  roots  that  represent  the  components  of  the  principal  stress  vector. 
The  components  are: 

0_ 
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-  •HI.41.  K;  (58) 

83  -  -2.41u  ^  • 

The  largest  difference  between  any  two  principal  stress  components  Is: 

I  ®l"®3l  ”  ^  •  (59) 

The  difference  between  the  radial  and  circumferential  stresses  Is: 


l'^rr“<^9el  "  ^  • 


(60) 


Accepting  Che  approximate  stress  state  (49)  and  estimating  the  shear 


strain  rate  Ou  /dr)  from  experimental  measurement  OU  /8r)  yields  Che 

Z  2 


observation  that  expressions  (59)  and  (60)  are  not  significantly  different. 

They  differ  on  Che  order  of  ten  percent  or  less  depending  upon  Che  relative 

magnitude  of  the  shear  strain  rate.  In  Region  III,  the  principal  and  general 

stress  states  are  Identical  by  definition  (dU  /dr  >  du  /dr  -  0).  In  Region 

z  z 

II,  however,  the  axial  deformation  Is  nonzero  but  as  just  mentioned,  Che 
distinction  between  the  expressions  (59)  and  (60)  Is  perhaps  not  significant 
for  practical  purposes.  Radial  flow  la  likely  to  predominate  In  both  regions. 
This  statement  Is  once  again  supported  by  experimental  and  nimerical  evidence 
(Glass  [64] ,  Pond  [69],  Wagner  [104],  etc.).  It  should  also  be  noted  Chat  Che 
yield  strength,  presumed  In  the  yield  criterion,  of  target  materials  subject 
Co  ballistic  Impact  Is  only  crudely  known.  The  errors  associated  with  the 
estimation  of  a  precise  value  most  likely  exceed  any  error  Introduced  by  the 
equlvalenclng  of  the  principal  and  general  stress  vectors.  With  this  In  mind, 
Che  stress  difference  term  In  equation  (44)  will  be  approximated  with  a 
principal  stress  difference  and  subsequently  Che  Tresca  yield  criterion  and 
the  Hohenemser-Prager  flow  rule  will  be  applied  in  the  same  manner  as  In 
Region  III.  The  resulting  radial  and  circumferential  stress  distributions  In 
Region  II  are: 
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for  a  <  r  <  a' . 


The  static  contributions  of  Che  general  stress  tensor  are  assumed 
negligible  since  Che  deformation  rates  and  viscous  stress  components  are 
estimated  as  large.  The  assumption  of  parallel  flow  eliminate!^  the  axial 
gradients  of  the  velocity.  The  axial  stress  component  Is,  therefore. 


negligible 
to  Che  axial  velocity  distribution  by: 


The  shear  stress  component, v  ,  Is  however,  finite  and  Is  related 
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where  the  gradient  term  may  have  Che  same  form  as  Che  experimentally  observed 
shear  strain  distribution.  In  order  Co  more  precisely  determine  Che  axial 
velocity  distribution,  equation  (45)  Is  recalled  and  Che  axial  gradient  of  Che 
axial  stress  is  Ignored: 
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Equation  (54)  crudely  estimates  Che  radial  gradient  of  Che  axial  velocity 
distribution  to  be  of  exponential  form,  with  a  decay  coefficient,  and 

maximum  shear  rate,  The  radial  derivative  of  this  distribution  will  also 

be  of  Che  same  exponential  form.  It  Is  assumed  that  Che  last  term  In  equation 
(64)  will  dominate  Che  right-hand  side  of  Che  equation  and,  therefore, 
equation  (64)  is  approximated  as: 


Providing  that  0(t)  Is  Independently  determined,  equation  (65)  is  separable 
and  can  be  solved  with  standard  methods.  Its  solution  is  given  in  Appendix  C 
With  this  distribution  in  hand,  the  complete  but  approximate  stress  dis¬ 
tribution  is  specified  in  terms  of  the  undetermined  differentiable  function, 
0(t). 

E.  The  Rod  Inversion  ReKlon  I 

The  steady  penetration  theory,  as  described  in  Appendix  A,  presumed  that 
the  penetrator/ target  interaction  surface  moves  along  the  symmetry  axis  with 
constant  velocity,  U.  If  a  coordinate  frame  of  reference  is  attached  to  the 
intersection  point  of  this  surface  with  the  symmetry  axis  (stagnation  point), 
then  penetrator  material  approaches  this  point  with  a  relative  velocity, 

S  -  V  -  U,  where  V  is  the  penetrator  velocity  in  the  laboratory  frame  and 
target  material  approaches  from  Che  opposite  direction  with  a  relative 
velocity,  U.  In  this  translating  coordinate  frame,  Che  axial  penetration 
process  appears  steady.  The  interaction  surface  has  a  characteristic  mushroom 
shape  Chat  remains  constant  throughout  the  steady  penetration  process.  Figure 
13  is  a  schematic  description  of  this  surface  along  with  a  hypothesized 
control  volume,  c. 
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Figure  13.  The  Penetrator  Interaction  Surface  and 
Control  Volume,  c 


The  penetrator,  characterized  by  its  circular  cross  section,  A^^ ,  density, 

P^,  and  velocity,  V,  upon  reaching  Che  stagnation  region,  inverts  and  flows 

along  Che  expanding  target  cavity,  through  the  incremental  control  volume,  c. 
As  in  Che  plastic  deformation  regions  of  Che  target,  incompressibility  is 
assumed,  material  enters  Che  control  volume  through  surface,  A^,  and  exits  Che 

control  volume  through  A^*  Assuming  that  Che  flow  is  steady  (no  acceleration) 

during  the  angular  deflection,  d0,  Che  target  must  supply  a  forced, 
determined  by  the  relation: 


which  Is  a  statement  of  the  conservation  of  linear  momentum.  The  inner 
surface  may  be  interpreted  as  a  free  stream  surface,  upon  which  the  pressure 
is  zero  and  from  Bernoulli's  equation,  the  speed  l?l  must  be  constant.  If 
mass  flows  through  surfaces,  and  A^  only  and  the  flow  velocity,  S  is 

uniform  over  each  of  these,  then  the  incompressible  continuity  equation 
requires  that: 

Ppl^l^l  “  fpl^|A2  *  Pp|s'lA3  »  constant  (67 


and,  therefore,  equation  (66)  becomes: 


F  *  Opl^l  sin  d0 


(68: 


Cross  section,  A^ ,  is  given  by  wr^  ,  where  r^  is  the  penetrator  radius.  The 


pressure  within  this  Incremental  control  volume  may  be  estimated  as  the  force 
acting  upon  the  outer  surface  over  the  area  of  this  outer  surface: 
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(69 


where  R  is  the  radius  of  curvature  of  the  outer  surface,  d0  is  the  angular 
deviation  of  the  velocity  vector  S  and  r  is  the  radial  position  of  the 
midpoint  of  the  outer  surface.  Substituting  equation  (68)  into  (69)  yields: 


1^  ^pITI  ^  sin  d0 
”  2  rRd0 


where  for  small  deflection  angles,  d0  >  sin  d0  and  then 
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The  steady  penetration  theory  requires  stress  continuity  at  the  stagnation 
point,  thereby  relating  IS<  and  U: 

Pp|3|2  -  pu2  (72 


where  p  is  the  target  density, 
written: 

P 


Equation  (7l)  may  then  be  alternatively 
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From  differential  geometry,  the  radius  of  curvature  is  related  to  the 
Instantaneous  surface  slope  and  its  derivative: 


(73 


d^z 

dr2 


In  order  to  evaluate  the  terms  dz/dr  and  d  z/dr  ,  the  steady  penetration 
theory  is  once  again  recalled.  The  time  rate  of  change  of  the  axial  position 
of  the  stagnation  point  Is  defined  as: 


dz 
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If  z>z(t(r))  where  t(r)  has  a  unique  Inverse.  r(t)  and  dr/dt  ^  0,  then: 

dz  _  dz  dt  _  „  dt 
dr  dt  dr  dr 
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dr  \  dr/  jj.2 


The  Inverse  functions  are  given  by  the  rules  of  differentiation. 


ii./i£l 

dr  \dt  f 


(H)^ 

and  when  combined  with  equation  (76)  and  written  In  more  compact  form  yield: 


Substituting  expressions  (78)  Into  equation  (73)  then  couples  the  pressure  In 
the  Incremental  control  volume  with  the  kinematic  properties  of  Che  outer 
surface  (r,  r,  r)  at  any  particular  time,  upon  any  axial  plane: 


prp2u3 


In  ocher  words.  If  the  radial  acceleration  process  of  the  target  cavity 
relative  Co  a  fixed  axial  plane  Is  known  a  priori,  Chen  the  time  variation  of 
the  pressure  In  Che  penetrator  at  the  axial  coordinate  can  be  determined. 
However,  It  Is  this  acceleration  process  Chat  is  the  desired  solution,  thereby 
suggesting  the  coupling  between  Regions  1  and  II. 


F.  The  Boundary  Coupling  of  Regions  I-IV  Yielding  the  Governing 

Differential  Equations 

Aa  noted  In  Section  IIA,  the  problem  is  divided  into  separate  regions  of 
different  deformation  behaviors  to  allow  an  analytic  formulation.  The 
formulation  of  Region  I  produces  a  relationship  between  the  pressure  within 
the  reversed  flow  penetrator  and  flow  properties  of  a  point  of  Intersection  of 
the  penetrator  Interaction  surface  and  a  particular  axial  plane.  Regions  II, 
III  and  IV  yield  the  stress  distributions  on  constant  z  planes  as  a  function 
of  radial  position  and  the  function  (l(t).  Regions  III  and  IV  presume 
negligible  axial  deformations  while  Region  II,  obeying  the  parallel  flow 
constraint,  allows  for  axial  velocity,  axial  displacement  and  shear  stress  in 
the  r-z  plane.  Before  the  stress  distributions  can  be  quantitatively 
described,  the  function  (9(t)  must  be  determined.  Indeed,  once  this  function 
Is  known,  then  the  displacement  and  velocity  fields  as  well  as  the  stress 
fields  will  be  determined. 

In  order  to  determine  the  function  0(t),  the  radial  stress  equations  In 
Regions  II,  III  and  IV  will  be  utilized.  The  boundaries,  whether  real  or 
Imaginary,  separate  these  Regions.  From  the  continuity  equation  applicable 
for  the  plastic  flow  Regions,  the  relation  for  the  radial  velocity  Is: 


If  the  velocity  of  a  particle  at  position  r,  at  time  t  Is  expressed  as  dr/dt, 
then  upon  Integration  of  equation  (80),  the  position  of  a  particle  originally 
located  at  r^  Is  given  by: 

r  -  [ro2  +  20(t)]^/^  (81) 

for  any  time  t.  This  expression  may  be  twice  differentiated  to  obtain  the 
acceleration  at  a  given  radial  position  and  time: 
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The  positions  of  the  boundaries  may  be  similarly  described: 


a(t)  -  [8^2  +  20(t)]^^^ 
a'(t)  -  [30-2  +  20(t)]^'^^ 

b(t)  -  [bo2  +  20(t)]^^^. 


(83) 


The  outer  boundary,  c(t)  can  be  considered  constant  If  the  target  Is 
assumed  seml-inf Inlte,  i.e.: 


c(t)  -  [cq2  +  20(t)3^^^ 
but  I Co2|  »  I 0| ,  so 


(84) 


C(t)  ^  Cq 


If  Che  lateral  dimension  Is  finite  and  the  outer  boundary  Is  likely  to 
permanently  deform,  Chen  the  elastic  Region  IV  vanishes  and  a  traction  free 
boundary  condition  can  be  applied  directly  Co  b(t). 


The  magnitude  of  the  penetration  cate  determines  the  nature  of  the 
matching  conditions  at  these  boundaries.  As  noted  in  the  introduction, 
framing  camera  experiments  involving  optically  transparent  materials  suggest 
that  shock  waves  are  likely  if  Che  penetration  rate  equals  or  exceeds  Che 
local  sound  speed  of  the  target.  If  the  penetration  is  "subsonic,”  shock 
waves  are  less  likely  and  Che  stress  and  velocity  distributions  are  then 
continuous  across  the  boundaries.  This  has  been  both  experimentally  and 
numerically  confirmed.  Pritchard  [70]  considered  subsonic  conditions  and  made 
measurements  at  gage  locations  on  the  penetration  trajectory.  He  found  chat 
as  Che  penetrator  interface  approached  the  gage,  the  stress  and  strain  field 
Increased  rapidly  (but  continuously)  from  zero  to  Che  stagnation  value  without 
any  evidence  of  shock  discontinuities.  The  computational  work  of  Wagner  [104] 
recognizes  the  same  behavior,  even  during  the  Initial  transient  phase.  For 
most  conventional  armor  applications,  the  sound  speed  in  metallic  target 
materials  is  typically  greater  chan  the  penetration  rates  of  both  kinetic 
energy  penetrators  and  shaped-charge  Jets.  For  these  situations,  stress  as 
well  as  velocity  continuity  will  be  enforced. 

The  radial  stress  variation  la  identical  in  both  Regions  II  and  III. 

Stress  continuity  at  boundary,  a' ,  Chen  allows  these  two  Regions  to  be 
considered  as  one.  This  combined  plastic  Region  Chen  has  an  inner  boundary, 
a,  and  an  elastlc/plastlc  outer  boundary,  b  (for  the  semi-infinite  problem). 
The  radial  stress  at  any  radial  position  r,  a£r<b  is  then  given  by  equation 
(61)  of  Region  II;  ”* 

a„<r)  -  arr(a)  “  {pS  +  ®y}  f  +  I  {p^^*  (85) 

The  condition  of  stress  continuity  on  the  Inner  boundary  Is: 

arr(a)  -  -P(t)  (86) 

where  P(c)  Is  determined  from  the  control  volume  analysis  of  Region  I. 
Combining  (86),  (85),  (83)  and  equation  (79)  of  Region  I  yields: 
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0(c)  Is  Che  undecermined  funcClon  of  Cime,  r  and  a  Che  inicial  radial 
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poslclons  of  Che  polnc  r(c)  and  Che  cavicy  wall  respeccively.  If  Che  scress, 
is  known  aC  a  polnC  r  for  all  Clme,  Chen  equaClon  (87)  becomes  a  second 

order  nonlinear  ordinary  differencial  equaclon  In  0(c).  Three  cases  In  which 
Chls  scress  boundary  condiClon  can  be  deCermined  are  now  considered. 

The  flrsc  case  involves  peneCraclon  InCo  mecallic  cylinders  of  small 
radial  dimension.  The  peneCraClon  process  Is  subsonic;  scress  and  veloclcy 
conClnulCy  Is  assumed  Chroughouc  and  Che  laceral  CargeC  boundary  Is 
uncons C rained.  This  siCuaclon  suggescs  boCh  Che  neglecC  of  Che  elasclc  region 
alcogecher  and  Che  Cracclon-f ree  condiClon  on  Che  ouCer  boundary.  Given 
,r^^(b)  *  0,  for  all  Clme,  equaclon  (87)  Chen  becomes: 
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The  second  case  Involves  Che  peneCraclon  Inco  very  massive  meCalllc 
Cargecs.  The  peneCraclon  process  Is  subsonic;  scress  and  veloclcy  conclnulcy 
Is  assumed,  buc  Che  CargeC  Is  seml-lnflnlce  or  of  very  large  laceral 
dimension.  An  elasclc  region  bounds  Che  plasclc  deformaclon  region.  Scress 
conclnulcy  ac  Che  elasClc/plasClc  boundary,  b(C)  requires  chaC: 

(for  all  Clme).  (89) 

Recalling  Che  Hohenemser-Prager  flow  condiClon  (equaclon  (26)  Secclon 
III),  applied  Co  Che  elasclc  boundary  polnC,  b. 
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and  remembering  ChaC  In  Che  elasClc  region,  (r)  >  -v  (r),  Chen 
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Equaclon  (87)  Chen  becomes: 
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The  Chlrd  and  final  case  Involves  deep  peneCraclon  Inco  low  sound  speed 
macerials.  The  peneCraclon  Is  supersonic;  a  shock  wave  represenclng  a 
dlsconClnulCy  In  sCress  and  parclcle  veloclcy  separaCes  qulescenc  maCerlal 
from  plasclcally  flowing  macerial.  If  Che  wave  is  assumed  Co  be  a  radially 


diverging  cylindrical  wave  Chat  follows  the  experimentally  observed  linear 
relationship  between  particle  and  wave  velocities,  then  the  radial  stress  may 
be  estimated  from  the  one-dimensional  Ranklne-Hugonlot  relations: 
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where  the  stress  Is  proportional  to  the  wave  speed  C,  the  particle  velocity  u 
and  the  undeformed  density,  p.  The  aforementioned  relationship  between  wave  ^ 
and  particle  velocities  is  given  by: 


Co  + 


Xu. 


(94) 


where  the  material  constants  are  available  in  tabular  form  from  several  shock 
property  compendia  [112,  113,  114].  Combining  equations  (93)  and  (94)  with 
the  particle  velocity  expression  (80)  at  shock  position,  b(t),  yields: 
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and  upon  combination  with  equation  (87)  the  governing  ordinary  differential 
equation  becomes: 
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6.  The  Determination  of  the  Initial  Conditions 


The  previous  section  coupled  the  deformation  Regions  by  applying 

appropriate  matching  conditions  at  the  boundaries  of  the  plastic  Regions  II 

and  III.  The  resulting  ordinary  differential  equations,  while  highly 

nonlinear  are  certainly  simpler  to  solve  than  the  original  non-linear 

partial  differential  equations.  Equations  (88),  (92)  and  (96)  are  not  likely 

to  yield  exact  solutions  In  their  full  form.  Either  additional  approximations 

must  be  applied  to  reduce  these  to  more  tractable  analytic  form  or  the 

equations  could  be  approximately  solved  with  many  standard  numerical 

techniques.  However,  before  either  approach  may  be  pursued,  two  Initial 

conditions  on  0(t)  and  the  Initial  values  a  and  b  must  be  independently 

o  o 

specified  to  make  the  system  determinant.  The  target  density,  P,  the 
penetrator  dimension,  r^,  the  steady  penetration  rate,  U,  the  target  yield 

strength,  the  target  viscosity,  ft,  and  the  shock  properties  C^  and  X  are 

problem  parameters  which  must  also  be  known  before  a  solution  Is  obtained. 
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The  cavity  expaoslon  process  begins  when  the  penetrator/ target  Interface 
arrives  at  a  particular  z-plane  of  Interest.  The  penetrator  Interface 
proceeds  through  the  remainder  of  the  target  In  a  steady  fashion  but  the 
trailing  cavity  expands  radially  through  an  unsteady  process.  It  Is  this 
acceleration  of  the  cavity  wall  In  the  radial  direction  that  Is  the  heart  of 
this  report.  Both  experiments  and  higher-order  numerical  solutions  have  shown 
that  there  Is  a  region  which  preceeds  this  moving  interface  that  accelerates 
the  originally  quiescent  target  material  to  some  undetermined  velocity  and 
stress  configuration.  Figure  14  shows  the  essential  characteristics  necessary 
to  establish  the  Initial  conditions. 


The  first  Initial  condition  has  already  been  Implicitly  defined.  The 
current  position  of  an  element,  a(t),  displaced  from  its  original  position, 
a^,  was  determined  from  the  Incompressible  continuity  equation  as: 

a(t)  -  [ao2  +  2B{t)]^^^.  ^ 

By  definition,  therefore: 


0(0)  •  0  . 
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The  Initial  radial  velocity  of  the  cavity  wall  at  this  Initial  radial  position 
has  also  been  determined: 

(99) 

but  at  this  point.  Its  value  Is  undetermined.  It  Is  known  that  the  value 
should  be  finite  and  positive.  Expression  (99)  along  with  the  requirement 
that  the  Initial  radial  velocity  be  finite  precludes  the  possibility  of  a^ 

equaling  zero.  Indeed  the  stress  distributions  trould  also  be  singular  upon 
the  axis  of  symmetry  (p>0)  if  this  were  allowed,  a  is,  therefore,  slightly 
offset  from  this  axis. 


Hill  [79]  and  Elchelberger  161]  recognized  the  necessity  of  determining 
the  radial  velocity  behavior  so  that  the  penetration  process  might  be  more 
completely  understood.  Hill  assumed  that  the  cavity  wall  Initially  maintains 
contact  with  a  symmetrical  ovoid  shaped  penetrator  ‘‘bulb'*  and  that  the  radial 
velocity  is  proportional  to  the  steady  penetration  velocity  and  the 
Instantaneous  slope  of  the  Interaction  surface  described  by: 


a  ■  A  sin  — 


The  radial  velocity  Is  then  defined: 
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*-«dI 


(100) 


(101) 


•  ir  „  A  /ir  z\ 
a  -  ^  U  -  cos  (3 


(102) 


For  the  initial  value,  z«0  and, therefore: 


4(0)  .  f  ^  u 


(103) 


but  the  ratio  A/L  Is  not  necessarily  known  a  priori.  If  the  surface  1: 
hemispherical  In  shape,  the  ratio  A/L  is  unity. 


Elchelberger  [61],  In  an  attempt  to  verify  the  proportionality  between 
a(0)  and  U,  performed  Instrumented  jet  penetration  experiments.  Considering  a 
steel  Jet  penetrating  a  steel  target,  a  Jet  collection  technique  together  with 
a  rotating  mirror  camera  and  a  film  switch/oscillograph  method  was  used  to 
simultaneously  detenalne  the  Jet  velocity  (V),  penetration  velocity  (U)  and 
mass  (M)  as  a  function  of  time.  The  incremental  mechanical  work  (dK)  done  by 
the  Jet  In  displacing  an  Incremental  volume,  dr,  was  derived  as: 


25  -  uv  ^ 

dr  dr 


(104) 


The  experimental  values  of  the  above  expression  were  plotted  against  U  and 
the  result  Is  presented  in  Figure  IS.  In  an  earlier  theoretical  effort,  Bethe 
[19]  derived  an  expression  for  the  work  required  to  expand  a  cylindrical  hole 
to  a  volume,  r.  in  differential  form  this  relationship  was  given: 


U  >  (cm/ jf  0} 


Figure  15.  Mechanical  Work  pec  Unit  Volume  Plotted  as  a  Function  of 
the  Square  of  the  Penetration  Velocity,  Eichelberger  [61] 
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where  ^  represents  an  average  flow  stress  during  plastic  deformation  and  b  was 

(  k\  ^ 

given  as  .455p  where  a  Is  the  radial  velocity.  By  measuring  the  slopes 

from  Figure  IS  directly,  a(0)  (the  Initial  radial  velocity)  can  be  estimated. 

For  penetration  velocities  greater  than  2.5x10^  cm/s,  the  radial  velocity  is 

given  by: 


a(0)  -  1.4U 


(106) 


and  for  penetration  velocities  below  2.5x10 


5 


It  appears  that: 


a(0)  -  .43U  (107) 

Is  more  appropriate.  Oplk  [39]  developed  an  approximate  hydrodynamic  theory 
for  the  motion  of  cavities  produced  by  the  Impact  of  meteorites.  By  assuming 
that  resistance  In  the  axial  and  radial  directions  Is  of  the  same  form,  he 
applied  the  steady  Bernoulli  equation  and  obtained  the  following  value  for  the 
cavity's  Initial  radial  velocity: 


a(0) 


(108) 


where  Is  the  penetrator  density, 
materials: 


a(0) 


For  similar  target  and  penetrator 


2 


(109) 


Considering  the  level  of  approximations  Invoked  to  obtain  the  above 
expressions,  the  precise  value  of  the  Initial  radial  velocity  Is  still 
unavailable.  To  at  least  narrow  the  limits  within  which  a(0)  may  vary,  an 
analogy  can  be  made  to  the  theory  of  compressible  axlsymmetrlc  stagnation 
flow.  White  [IIS]  discusses  the  case  of  fluid  flow  past  bodies  of  revolution 
of  finite  dimension,  D,  and  approach  velocity,  U.  The  value  of  the  local 
velocity  gradient  K  at  the  stagnation  point  and  the  distance,  x,  from  this 
point  determine  the  velocity  of  the  flow  field  on  the  edge  of  the  boundary, 
Og.  a»: 


The  velocity  gradient  K  depends  upon  D  and  U  and  to  a  lesser  extent,  upon  the 
body  shape  and  approach  Mach  number.  For  subsonic  fluid  flow  (Mach  number 
<1),  the  Raylelgh-Janzen  procedure  can  be  used  to  estimate  the  velocity 
gradient  R.  Two  body  shapes,  presented  by  White,  are  suitable  for  further 
Investigation.  They  are  the  flatnosed  cylinder  and  the  sphere.  Experiments 
conducted  by  Fugelso  [63]  and  Welhrauch  [73]  consider  the  penetration  of  high 
speed  metallic  rods  (V  >  l'-2  km/s)  Into  thick  finite  targets.  The  shapes  of 
the  interaction  surface  are  observed  with  flash  x-ray  photography.  The 
surfaces  nearest  to  the  stagnation  region  appear  to  be  somewhat  between  flat 
and  hemispherical.  At  higher  velocities,  shaped  charge  Jets  have  been 
radiographed  (unpublished  Ballistic  Research  Laboratory  experiments)  while 
penetrating  into  homogeneous  targets.  The  penetrator  Interaction  surfaces 
have  similar  shapes  as  those  of  the  constant  velocity  rod  experiments.  Finite 
difference  calculations  of  Harlow  [55]  and  Kucher  [56]  further  support  these 
observations. 
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Throughout  this  report,  It  has  been  assumed  that  the  penetration  process 
is  steady  when  considering  the  axial  penetration  rate.  The  accuracy  has  been 
^perlmentally  and  numerically  confirmed.  In  addition  to  the  rate  being 
steady,  the  shape  of  the  Interaction  surface  is  assumed  to  remain  unchanging 
as  the  penetration  proceeds.  The  validity  of  this  approximation  is 
qualitatively  confirmed  by  the  experiments  just  mentioned.  If  the  target 
material  behaves  like  a  fluid  in  this  stagnation  region,  then  the  expressions 
for  the  sphere  and  blunt-nosed  cylinder,  presented  by  Vhlte  [115],  determine 
the  initial  radial  velocity  of  the  cavity  wall  near  the  axis  of  synmetry. 
Recalling  equation  (110)  the  radial  velocity,  a(0)  is  related  to  the 
stagnation  gradient,  K  and  the  initial  radial  position,  a^  by: 


If  the  initial  radial  velocity  a(0)  occurs  at  the  point  where  surface 
curvature  begins  (up  to  this  point  the  surface  is  flat)  then  the  blunt-nosed 
cylinder  configuration.  Figure  16,  is  plausible.  White  notes  that  while  the 
compressible  stagnation  gradient,  K  is  not  known  for  this  configuration,  it 
should  be  less  than  the  incompressible  spherical  and  greater  than  the 


supersonic  blunt-nose  measurement  IK  >  .15 —  I  of  Trimmer  [116].  White 

\  *o/ 

U 

determines  an  incompressible  value  of  vlk  —  ,  equation  (111)  then  becomes: 


Figure  16.  Blunt-nosed  Cylinder  Inan  Uniform  Flow 


Tf  the  interaction  surface  geometry  la  more  hemispherical  near  the  axis  of 
an"lo:  is  subsonfc  (0<.6  C„).  than  the  Raylel*h-.e„ten 

procedure  eay  be  applied  to  eatlute  the 

fu71  cooalder.  •  velocity 

solution  Is  expanded  in  powers  .  r  _  -  snhere  of  diameter,  D  and  flow 

and  pressure  distributions  are  presented  for  a  sphere  oi  niamei.  , 

velocity,  Ue 


tlista  (d««r««s) 


thata  (daflraas) 


Figure  17.  The  Velocity  and  Pressure  Distributions  upon  the 
Surface  of  a  Sphere  (Kaplan  lll7j) 

Near  the  stagnation  point  (fl-0)  the  velocity  varies  linearly  with  radial 
CdS6  * 

-  3  (1  -  0.252Mo2  -  0.175Mo^  +  .  .  .)  (li; 


where  the .Mach  number  Is  defined  as  Che  ratio  of  the  velocity  U  to  the 

characteristic  wave  speed,  C^.  From  radiographic  evidence  it  can  be  observed 

Chat  the  diameter  of  Che  Interaction  region  Is  typically  less  Chan  four  times 
Che  penetrator  radius,  r^,  or 


D  <  4  Tp 


and  upon  combination  with  equations  (111)  and  (113)  yields, 
a(0)  >  T  —  U  (l  -  0.252Mo2  -  0.175Mo^  +  .  .  . 


(114) 


(115) 


This  inequality  is  valid  only  for  penetration  rates  less  chan  .6  C^.  As 

an  example  of  its  utility,  however,  consider  Che  case  of  a  kinetic  energy 
penetrator  perforating  a  steel  target  at  2  km/s  (V  «  3.5  km/s).  The  sound 
speed  is  approximately  6  km/s  and  Che  Mach  number  is,  therefore,  .33.  The 
difference  between  Che  Incompressible  and  compressible  estimates  is  only  three 
percent.  This  difference  will  become  greater  as  the  Impact  velocities  increase. 
The  ratio  undetermined  but  should  be  of  Che  order  of  unity 

if  the  limit  (115)  is  to  agree  reasonably  with  the  earlier  estimates  of  Oplk, 
Elchelberger  and  White  (blunt-nosed  cylinder).  The  investigations  above 
suggest  the  following  approximation: 

ao  "  (116) 

where  the  constant  K  Is  between  0.4  and  ‘JT  for  geometries  somewhere  between 
flatnosed  and  hemispherical  and  for  all  penetration  velocities  (compressible 
or  not).  If  the  constant  K  can  be  presumed  to  be  known,  then  two  Initial 
parameters,  a^  and  b^,  still  remain  to  be  determined  before  the  cavity  growth 

equation  (87)  may  be  integrated. 

The  two  initial  conditions,  (}(0)  and  (9(0)  are  restated: 


0(0)  -  0 


(117) 


and  from  equations  (99)  and  (116): 


0(0)  -  aoKU 


(118) 


The  maximum  pressure,  P^,  occurs  at  the  moving  stagnation  point  along  Che 

symmetry  axis.  Kaplan  [117],  considering  subsonic  flow  around  a  sphere, 
observes  that  Che  pressure  varies  only  slightly  from  the  stagnation  value  for 
angles  up  Co  approximately  ten  degrees  from  Che  symmetry  axis.  Therefore,  for 
initial  cavity  positions  close  Co  the  symmetry  axis,  Che  pressure  is  given 
approximately  by  Che  Bernoulli  equation  (Appendix  A)  as; 


•|  P 


(119) 


where  Che  density  is  Che  ambient  target  density  and  U  represents  Che  steady 
penetration  velocity.  The  pressure  in  Che  penetration  Region  I  was  derived 


MO  .  i  .  ^ 


02/a| 

u2a2\3/2  * 


^  If) 


(120) 
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The  value  of  P  at  time,  t  >  0  Is  determined  by  substituting  (98)  and  (118) 
Into  expression  (120)  and  upon  equating  this  value  with  equation  (119),  the 
following  relation  results: 


3/2  ,.2 


liio  -  kVI  "  11  + 


(121 


where  9^  represents  the  Initial  value  of  (}(t).  The  absolute  value  character 


of  (121)  can  be  eliminated  by  the  following  argument.  The  expression  for  the 
radial  acceleration  of  a  point  on  the  cavity  surface  was  determined  as: 


...  ,  0(t)  0^(t) 

"  a(t)  a(t)3  * 


(122 


Its  Initial  value  at  the  radius  a  ,  with  the  assumed  Initial  conditions  is: 

o 


ao  ”  ~  1^0  “  , 


(123: 


From  inspection  of  the  shape  of  the  interaction  region  during  steady 


penetration,  the  curvature  requires  that  a^  be  negative.  Since  a^  is  finite 


and  positive,  equations  (123)  and  (121)  then  require  that 

2 


00  -  !  k2  -  [1  +  ^  U2  . 

'■p 


(124: 


The  radial  equation  of  motion  along  with  the  condition  of  stress 
continuity  at  the  cavity,  at  time  t  •  0  is: 


p  u' 


(125: 


for  the  subsonic,  finite  lateral  dimension  case. 


P0O  —2  i 


bo  2 


(126 


for  the  subsonic,  seml'lnf Inlte  lateral  dimension  problem  and 

2 


P0O  +  Oy[  4n  ^  +  \\o  k2u2  -  (l  -  P  ^  KU  (c^  +  X  ^  KU 


b  2, 


2  ^ 


(127) 


for  the  case  of  supersonic  penetration. 

Combining  equations  (124),  (125)  and  (126)  dividing  through  both  sides  by 


2  2 
pU  and  neglecting  the  small  terms,  v  /pU  and  (a  /b  )  yields: 


which  Is  approximately  true  for  all  three  cases.  If  further, 

u/KpaoU  «  1 


(128) 


(129) 


an  assumption  which  may  be  verified  by  applying  the  material  constants  p,  p, 
the  penetration  velocity,  U  and  recognizing  that  a  “O  (r  ),  then  equation 
(128)  reduces  to:  o  P 


k2  -  [1  +  ^  _  i  +  K2y 


(130) 


For  the  first  boundary  case  considered,  the  outer  boundary,  b^,  is  known 

and  equation  (130)  relates  a^  with  the  estimated  value  of  K.  For  the  other 

cases,  however,  the  boundary  location,  b^,  is  not  necessarily  known  a  priori. 

In  order  to  determine  a  and  b  for  these  two  cases,  a  control  volume  is 

proposed,  centered  on  the  radial  position,  b^,  moving  with  the  boundary 

velocity  b  .  Ignoring  the  negligible  curvature  as  the  radius  increases,  the 
one  dimensional  conservation  of  radial  momentum  requires  that: 


jj.(bo)A  -/*bo  P  bods  -  pbo^A 


(b„)  -  p  ?  -pkV  . 


vr''°0'' 


(131) 


Equating  (131)  with  the  stress  on  the  elastic/plastic  Interface, 


aK2u2  -  I  (oy  +  u  \  ^ 


2  2  ry 


(132) 


and  once  again  recognizing  that  a^/b^«l  then. 


fo  ^  (av/2Pu2y^^ 


(133) 


Combining  equations  (133)  and  (130)  establishes  a  relationship  between  the 
initial  cavity  radius  and  the  initial  cavity  velocity  coefficient,  K: 


k2  -  [1  +  k2]^/2  ^  I  In  ^  _  i  (1  +  k2)  . 

*  I'd  I 


(134) 
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Given  a  value  of  K,  the  Initial  radius,  a^,  the  initial  plastic  boundary, 

and  the  Initial  condition  0(0)  can  be  estimated.  This  estimate  may  be 
Improved  by  relaxing  the  approximation  (129)  and  iteratively  adjusting  a 
until  equation  (128)  is  satisfied  to  the  desired  accuracy. 


III.  SOLUTION  OF  THE  INITIAL  VALUE  PROBLEM 

The  preceding  chapter  established  an  iterative  method  for  determining  the 
initial  conditions  necessary  for  integration  of  the  governing  equations; 
either  exactly  or  numerically.  The  method  required  that  the  Initial  radial 
stress  equal  the  stagnation  pressure  at  an  Initial  radius  slightly  displaced 
from  the  symmetry  axis  and  that  the  radial  stress  at  the  plastic  boundary  be 
Initially  given  by  the  one  dimensional  momentum  jump  condition.  The  cavity 
then  radially  expands  in  a  decelerating  fashion  until  the  stresses  in  the 
plastic  regions  reach  values  of  the  order  of  the  static  yield  strength.  When 
this  happens,  the  problem  is  for  practical  purposes  completed.  It  is  quite 
likely,  however,  that  elastic  deformations  will  continue  until  the  stresses 
vanish  completely;  but  due  to  their  relative  magnitudes,  they  may  be 
practically  neglected.  The  problem  will  be  assumed  complete  when  the  stresses 
everywhere  in  the  plastic  regions  reach  the  static  elastic  range. 

A.  Runge-Kutta  Solutions  to  the  Complete  System 


The  integration  of  ordinary  differential  equations  by  numerical  techniques 
is  well  established.  Choice  of  the  specific  method  is  usually  based  upon  the 
ease  of  application  as  well  as  the  type  of  equation  considered.  In  general, 
equations  of  second  or  higher  order  can  be  reduced  to  systems  of  first  order 
equations,  e«g., 

y  >  f(x,  y,  y)  can  be  written  as 
y  -  z 

z  -  f(x,  y,  z) , 

a  pair  of  simultaneous  first  order  equations.  One  of  the  most  common  methods 
used  to  solve  this  system  of  first  order  equations  is  the  process  of  Gill 
[118].  It  differs  slightly  from  the  fourth  order  accurate  method  originally 
proposed  by  Runge  [119]  in  the  precise  values  of  the  expansion  parameters 
(Cohen  [120]  compares  the  two  methods).  A  FORTRAN  program  of  Gill's  method  is 
given  as  an  appendix  in  the  text  by  White  [115]  and  reproduced  here  (Tables  1 
and  2)  with  the  programming  logic  appropriate  to  this  problem. 

The  logic  located  between  statements  10  and  20  are  expressions  for  the 
system  of  first  order  differential  equations  representing  the  second  order 
differential  equation  in  0(t)  for  the  subsonic,  semi-infinite  target  case 
(92).  The  input  parameters  are  described  in  order  of  the  occurrence: 

X  >  initial  problem  time  (zero  for  this  problem) 

XLIM  >  maximum  problem  time  (limits  actual  computer  run  time) 

H  «  time  step  (constant  for  this  problem,  but  not  restricted) 

M  »  an  Index  which  must  be  set  to  zero  by  the  programmer 
N  ■  number  of  first  order  differential  equations  (2  for  this  program) 

2 

SIGY  >  Invlscld  target  yield  strength  (dynes/cm  ) 

C0  ■  bulk  sound  speed  of  the  target  (cm/s) 


L-n  1.-4. 


c 

c 

c 


Table  1 .  Runge-Kutta  Program 


PRQ«>AII  SCaTT(XNPUT*0Ur»UT>TAPE9>INPUT«TAPE6«DUTPUTI 

aTN  ORDER  RUNDE  KUTT*  SOLUTION  I  DILL'S  NETHOO  I 

OINEMSXON  Y(4|«R(4) 

REU  NU«LANOA 

1  E0RNAT«2X.*5XDT«»»E10**»*C0«*»EI0.3»»RH0T«*»FX0.l»'  NU»**E10.2» 'R 

tF«'*fX0.2«  •iaN0a***E10.2»*U»»»EX0.E/E*»*n«*#E10.Sf  •T2»*»Ei0.3»2EX 
10.21 

4  raRNAri3FX0.4«22S> 

3  F0RNATCRFX0.3) 

9  FORNATIIOX**!*  •*EX0.3»*XLIN-  •*EX0.3»>H«  '.EXO. 3*2X5) 

XI  FORMAT (2X* ‘TINE-  'tEXO.!** FHXOOT-  •*£X0.3»*FHX>  •»EX0.1#»4-  ‘.EXO. 

13*'FeHECX«  *«EX0.3» 

S3  FQRNATISEXO.S) 

REaOIS*4IX*XLXN«H*N*N 

REAO(S»SISXSY*CO«RHOT,IW>RF*R*U*LANOA 
REA0(3*S)Y<X)>Yf 2)>D0*A0 
MRXTE(4«9IX«XLXN*H*H*N 

IIRITEIA»3ISXSY>CO*RrtOT*ltU»RF*LANDA»U*Y(l),Y(2)*BO*AO 

JFLAD-0 

UO.O 

ZX-0.0 

a  XFII-XL1H)****7 

4  CALL  RUNDE<N*Y*F*X*HfN*KI 
DO  TO  (X0f20)«K 

XO  Ala A04«2.«2. 99(2) 

AaSORTCAXI 

DXaa09«2.92.9Y(2) 

DaSORTfDXI 

OaO. 5*  RHaT«RF«92 ••UFFS. 

CaUOD<A/DI 

A2a(i.aU9«2.«U/Y(Xt992.)99X.S 

TXaOF(A*Ya)«A2) 

Aaa<X.FIXI-(X.FAXI 

T2a<  {RH0T«Yai«92.  )<-ffU*Y(l  t  IRAS/Z. 

Taa0.9«<SIDYaHU*Y(XI/lll 

TAaSXDYK 

FXaF (1 )• A/ (Y(1I««3.*A2) 

F2al./tA9Ylll*A2) 

FCHECRaO^CFX-FZt 

IF(JFLaD.ED.OI  FCMECRaO.S9RHOT*U9*2. 

JFLADal 

FaAaStFCHECX) 

T0aRH0r*C 

F < X I a< T2aT 3-T4-F » FTO 
IF(Y(X)»LE.O.)  DO  TO  7 
FIEJaTIll 
DO  TO  4 

20  ZaS9RT(A0992.92.9YUn 
ZZaZZ^l. 

ZFRXNTalOOO. 

XFfZl.ED.O.I  DO  TO  2S 
XFIZZ.NE.ZFRINT)  DO  TO  8 
Z9  VRXTE(4*33)AX«A«8l*D*A2*A3*rO*rX 

VRXTE(4»33)TZ*T3»T4»F(1)>F(2I*F0»FX*F2 
VRITE(4*XX)X*Y(X)*Yf2l*Z»FCNECR 
ZZ-0.0 
zx-x.o 
DO  TO  a 
7  STOP 


RHOT  >  target  density  (g/cm  ) 

2 

MU  ■  target  viscosity  (dyne  -  s/cm  ) 

RP  -  penetrator  radius  (cm) 

U  »  steady  penetration  rate  (cm/s) 

LAMDA  >  slope  of  the  US-UP  data  (nondlmenslonal) 

Y(l)  ■  0(0)  Initial  condition  (cm^/s) 

Y(2)  -  0(0)  initial  condition  (cm^)  - 

B0  -  initial  position  of  plastic  boundary  (cm) 

A0  ■  initial  position  of  cavity  surface  (cm). 

Table  2.  Subroutine  Runge 


•XlliS  NfTH00-«UNfiE  KUTTA 


X 

2 


t 

4 

9 


SUMWTINE  RUN6EIN,r»P«|,H*ll,K| 
OtNENSZON 

M  TO  ll*4«5*3*7t>N 
00  2  t>l*N 
0(tl«0. 

AaO.S 
CO  TO  4 

A>1.70n04TCilC4 3479244 
2a|*H*0.9 
00  4  I>UN 

T(t|arai«A*<Eat*H-4ai| 

4  C(II«2«eA4H*E(I}^IX.>2»4A)4C(X> 
A««242  U21iO 194924794 
CO  TO  4 
T  00  C  I«1«N 

c  r(i}«r(xi«M»9ai/4.«c(X)/9. 

HpO 

K«2 

CO  TO  10 
4  R>1 
10  RCTUm 
INO 


The  case  of  subsonic  penetration  into  a  target  with  small  lateral  dimension 
can  be  considered  by  setting  the  program  variable,  T3,  to  zero  throughout  the 
calculation.  The  program  then  solves  equation  (88)  which  was  derived  by 
applying  a  traction  free  boundary  condition  to  this  lateral  surface.  The 
supersonic  penetration  case  can  be  similarly  considered  by  substituting  the 
expression  for  T3  with, 

T3  -  RHOT*(Y(1)/B)*(C0  +  LAMDA*Y(l)/B) 

which  represents  the  Jump  stress  at  the  shock  surface  (also  the 
elastic/plastic  boundary)  derived  from  the  one-dimensional  Rankine-Hugonlot 
relations. 

The  material  properties,  SIGY,  C0,  RHOT,  LAMDA  and  MU  must  be  supplied 
since  no  default  values  are  included.  The  first  three  may  be  found  in  con¬ 
ventional  material  handbooks.  The  viscosity  coefficient,  MU,  may  be  estimated 
in  many  cases  from  the  literature  review  of  Walters  [125],  as  described  in 
Appendix  B;  otherwise,  an  independent  experiment  is  required  to  determine  this 


parameter.  The  design  parameters,  RF  and  U,  must  also  be  specified  before  the 
calculations  may  proceed.  In  the  case  of  a  kinetic  energy  projectile,  RF  Is 
the  design  value,  but  for  high  velocity  jets,  a  precise  value  is  not  usually 
known.  Estimates  can  be  made  from  flash  x-ray  photographs  taken  of  the  jet 
before  and  after  It  Interacts  with  a  particular  target  element.  The  pene¬ 
tration  velocity,  U,  can  either  be  determined  experimentally  or  calculated 
from  the  one  dimensional  steady  model  (Appendix  A). 


The  Initial  conditions,  Y(l),  Y(2),  A0  and  B0  are  entered  as  input  and 
must  be  determined  by  the  method  described  in  the  preceding  chapter.  The 
remaining  Input  parameters,  X,  XLIM,  H,  M  and  N  are  related  to  the  solution 
execution.  The  Initial  time,  X,  Is  zero.  The  solution  proceeds  out  to  the 
time  specified,  XLIM,  or  to  the  point  where  the  velocities  drop  to  zero  (all 
stresses  are  of  elastic  magnitude) .  M  and  N  should  be  set  to  zero  and  two 
respectively.  The  time  step,  H,  Is  constant  throughout  the  calculation.  Its 
value  must  be  selected  with  care.  The  characteristic  time  of  most  cavity 
expansion  processes  la  of  the  order  of  microseconds,  consequently  H  must  be  of 
an  order  less  than  this  to  Insure  temporal  resolution  as  well  as  numerical 
convergence.  Since  the  length  of  the  calculation  depends  inversely  upon  the 
step  size,  an  Infinitesimally  small  value  Is  prohibitive.  A  suggestion  for 
the  determination  of  an  appropriate  step  size  Is  that  one  selects  a  reasonably 
-9 

small  value,  say  10  seconds  and  allows  the  program  to  run  Its  course.  Then 
select  an  order  of  magnitude  smaller  value  and  let  the  program  run  once  again. 
The  difference  between  the  two  calculations  should  be  very  small  and  the 
solutions  should  converge  as  the  step  size  Is  progressively  diminished, 
choice  of  the  final  step  size  will  then  be  limited  between  a  value  that 
ensures  numerical  convergence  and  a  value  that  is  reasonable  from  a  cost 
(computer  time)  viewpoint.  For  the  range  of  problems  considered  in  this 


The 
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report,  a  value  of  10  seconds  Is  appropriate. 


B.  Re4uctlon  to  Quadrature  via  Additional  Approximations 


The  equations  developed  In  Section  II  represent  the  cavity  growth 
behavior,  stress  distributions  and  velocity  distributions  that  can  be  expected 
during  deep  steady  penetration.  Approximations  concerning  the  nature  of  the 
behavior  were  made  after  reviewing  experimental  and  computational  evidence. 

The  governing  equations  of  motion  were  subsequently  reduced  to  single,  second 
order,  ordinary  differential  equations  which,  unfortunately,  are  highly 
nonlinear,  thereby  making  their  exact  solution  very  difficult.  If  not 
Impossible.  In  Section  IIIA  a  Runge-Kutta  technique  Is  described  which  can  be 
utilized  to  solve  the  full  equations  numerically.  In  this  section,  additional 
approximations  are  Imposed  to  simplify  the  equations  so  that  an  exact  solution 
Is  possible. 


Suppose  the  viscosity  Is  negligible  and  the  pressure  at  the  Inner  boundary 
vanishes  quickly  so  that  the  cavity  expansion  process  Involves  only  the 
Initial  momenta  and  stress  distributions,  then  equation  (92)  reduces  to: 


P0  + 


1 

-  2  Oy. 


(135) 


The  plastic  boundary,  b.  Is  typically  much  greater  than  the  cavity  dimension, 
a,  thereby  suggesting  the  approximation: 


If  further  It  Is  assumed  that  the  logarithmic  term  does  not  vary  significantly 
during  the  expansion  process. 


A  *0 

In  ^  ■  In  ^  ■  Ki 
b 


then  equation  (135)  slmpllflea  to: 


i  +  i  —il - il-i-.  i! 

2  Ki(ao2+20)  P  1 2Kl  » 


(137) 


(138) 


Define  the  new  variables, 

y  ■  "I  Ki(ao2  +  20)  “  n(ao2  +  20) 

y  ■  2n0 


y  ■  2n0 


■  p  \  n  y 
then  (138)  becomes: 


(^) 


2nAiy  -  yy  +  . 


(139) 


Ut  P  ■  J  and  y  -  ji  p  -  jS  P  .  P  g 


then  (139)  becomes: 


dP  p2 

2"Aly  -  yp  2^^  * 


(140) 


Let  2  ■  P  then  —  ■  2P  —  and  (140)  becomes: 

4nAi  “  ^  • 

*  dy  ny 


A  Ay  ■  I  w  2 

"  iiy  ^  “  n  - ^2 -  then  (141)  becomes: 


(141) 


j  .  .  I 

4nAi  ■  ny  +  nm  +  m 


or  rearranged  as: 


dm 


y  4nAi  -  m(n  +  1) 
Integrating  both  sides  yields: 


or 


yCi  ■  <4nAi  -  (n+l)m 


n 

n+1 


or 


4nAi  -  (yC]^) 


n+1 

n 


n+1 


(142) 
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(144) 


Performing  the  Inverse  transformations,  z  ••  my,  P  >  Vz  and  4^  >■  P  yields: 

^  dt 


dt 


•1/2 


dy. 


(1A5) 


The  problem  has  been,  therefore,  reduced  to  numerical  quadrature.  A  function, 
t  •  f(y)  will  result.  The  Inverse  of  this  expression,  y  ■  f(t)  then 
determines  0(t)  through  the  Inverse  transform: 


0(t)  - 


(146) 


The  advantage  of  this  approach  over  the  numerical  solution  of  the  full 
equations  Is  not  obvious.  Vfhlle  the  quadrature  may  be  performed  by  hand,  a 
more  accurate  and  timely  approach  would  be  to  utilize  a  computer.  The 
difference  In  computer  running  cost  between  a  Runge-Kutta  solution  and  a 
finely  partitioned  Integration  algorithm  Is  most  likely  negligible. 
Furthermore,  the  accuracy  of  the  quadrature  solution  Is  limited  by  the 
approximations  Invoked  to  reduce  the  governing  equations  to  the  form  (138); 
these  approximation  errors  typically  exceed  any  error  Introduced  by  the  fourth 
order  Runge-Kutta  algorithm. 

Alternate  sets  of  approximations  can  be  applied  to  the  full  equations 
(88),  (92)  and  (96)  reducing  them  to  varying  forms,  some  of  which  may  have 
analytic  solutions.  The  overwhelming  limitation  to  this  approach,  however,  is 
In  the  ability  Lo  estimate  the  errors  Introduced  by  each  approximation. 
Therefore,  for  the  remainder  of  this  report,  the  numerical  solution  to  the 
full  equations  will  be  pursued. 


IV.  PARAMETRIC  DEPENDENCE  AND  COMPARISON  WITH  EXPERIMENTS 


Section  II  dealt  with  the  formulation  of  the  cavity  expansion  problem. 

The  generally  complex  system  of  governing  field  equations  was  reduced  to 
Initial  value  problems  by  applying  approximations  suggested  by  experiments  and 
previous  numerical  Investigations.  Appropriate  Initial  conditions  were 
determined  from  problem  parameters  by  specification  of  the  boundary  stress 
conditions.  Section  III  dealt  with  possible  solutions  of  the  Initial  value 
problems.  A  Runge-Kutta  technique  was  proposed  to  solve  the  full  equations 
(88),  (92)  and  (96)  with  properly  specified  initial  conditions.  The  full 
equations  were  reduced  to  a  simpler  form  only  after  additional  approximations, 
consequently  the  accuracy  of  the  simplified  solution  was  in  question  and  the 
advantage  to  be  gained  by  the  simplification  was  reduced.  This  section  will 
present  Runge-Kutta  solutions  to  the  full  equations  as  developed  in  Section 
IIF. 


The  governing  equations  and  initial  conditions  Involve  the  previously 
defined  problem  parameters:  p,  »  Vy*  A,  C^,  U,  r^,  K,  a^  and  b^.  The  Impact 

velocity,  penetrator  density  and  penetrator  strength  characteristics  are 
Implicitly  Included  within  the  penetration  velocity,  U  (as  described  in 
Appendix  A).  The  functional  Influence  of  these  parameters  on  the  growth 
behavior  of  the  target  cavity  appears  explicitly  In  equations  (88),  (92)  and 
(96)  but  the  relative  Importance  of  each  parameter  for  the  prediction  of  the 
final  cavity  dimension  or  the  Instantaneous  growth  rate  has  not  yet  been 
determined.  Within  the  following  several  sections,  the  parametric 
dependencies  will  be  presented.  The  particular  examples  selected  are  those 
that  have  been  studied  experimentally,  thereby  allowing  a  comparison  between 
the  model's  predictions  and  actual  measurements.  Unfortunately,  the  growth 
rate  cannot  be  observed  experimentally.  However,  the  time-integrated  final 
value  can  and  will  be  utilized. 

A.  Dependence  of  the  Initial  Conditions  (K,  a^,  b^) 


As  already  noted  In  Section  IIG,  the  precise  value  of  the  initial  radial 
velocity,  a(0)  Is  not  known.  Estimates  are  available  from  different  models 
and  the  range  within  which  a  (0)  Is  expected  to  vary  Is  0.4U  <  a(0)  £  1.4U. 

If  K  Is  known,  then  a^  and  b^  are  determined  by  the  iterative  method  described 

In  Section  IIG.  Therefore,  K  will  be  chosen  between  these  limits,  a  and  b 

’  o  o 

will  be  subsequently  determined  and  the  growth  behavior  will  be  predicted  for 

the  following  situation. 

Silsby  [121]  and  Blanks  [122]  performed  experiments  that  involved  the 
normal  Impact  and  penetration  of  long  tungsten  alloy  rods  into  armor  steel 
targets.  Impact  velocities  were  varied  between  2  and  4.5  km/s  and  the  final 
hole  diameters  were  measured.  The  resulting  cavities  were  cylindrical  in 
shape  and  the  quoted  final  cavity  diameters  represented  averages  of  the 
diameter  at  several  depths.  The  lateral  dimensions  of  the  steel  plate  were 
quite 'large  and  the  sound  speed  In  steel  is  approximately  6  km/s,  therefore 
equation  (92),  the  subsonic,  semi-infinite  case  Is  selected  to  model  this 
problem.  The  following  values  are  assumed; 


p  >  17.3  g/ cm  (measured) 

P  3 

p  >  7.85  g/cm  (measured) 

r  «  0.39  cm  (measured) 

P  10  2 

<r  “lx  10  dynes/cm^  (GoldsmiCh  [3]) 

y  5 

V  -  4.4  X  10  cm/s  (two  measured  values 

V  -  3.3  X  10^  cm/s  are  selected  for  study) 
and  the  viscosity,  /i.  Is  given  by  Harlow  and  Pracht  [55]  as 

4 

/I  «  2  X  10  poise  (fit  to  the  data  of  Walsh,  et  al. 
[123]). 

The  penetration  velocity  Is  then  determined  by  equation  (A.3)  for  the  two 
chosen  Impact  velocities.  Representative  values  for  the  constant  K  are 
selected  and  corresponding  and  ate  determined.  Figure  18  presents 

the  Runge-Kutta  solutions  for  the  4.4  cm/s  Impact  velocity  and  the  following 
Initial  condition  parameters: 


selected  values 


1.400 


K  -  .700 


values  Iteratlvelv  determined 


a  •  .205  cm  b  <■  3.180  cm 
o  o 

a^  «  .132  cm  b  •  1.150  cm 
o  o 

a  *  .112  cm  b  •  .838  cm 
o  o 


Values  of  K  less  than  '».70  were  unable  to  satisfy  equation  (128)  for  any 
Initial  cavity  dimension.  Therefore,  the  correct  solution  should  lie 
somewhere  between  the  predictions  corresponding  to  K  -  1.4  ^md  K  *  .70.  The 
value,  .78,  represents  the  Incompressible  blunt  cylinder  configuration 
presented  by  White  [115].  The  predicted  final  cavity  radius  for  this  value 
agrees  quite  well  with  the  experimental  measurement  of  Sllsby  [121].  The 
final  hole  dimension  does  not  scale  linearly  with  the  constant  K,  although  the 
magnitude  Increases  as  R  Increases.  The  time  required  for  the  stresses  to 
decay  to  elastic  magnitude  Increases  with  Increasing  K. 

Figure  19  presents  the  Runge-Kutta  solutions  for  the  3.3  km/s  Impact 
velocity  with  the  following  Initial  condition  parameters: 


selected  values 


K  -  1.400 


K  -  .780 


values  Iteratively  determined 


a  *  .200  cm  b  >  2.220  cm 
o  o 

a  ~  .119  cm  b  >  .732  cm 
o  o 


Once  again  the  value  of  K  «  .78  yields  a  final  cavity  dimension  that  agrees 

reasonably  well  with  the  Independent  experiment  for  V  -  3.3  km/s.  Comparing 

Figure  18  and  19  for  the  same  K,  the  diameter  of  the' cavity  and  the  time 

required  to  reach  the  final  value  Increase  with  Increasing  Impact  velocity. 

Additional  comparisons  with  experiments  for  variable  Impact  velocities  will  be 

given  In  the  next  section.  In  addition  to  the  cavity  growth  curves,  the 

stress  distributions,  «  and  <r..,  are  presented  In  Figure  20  for  three  times 
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Parametric  Dependence  of  the  Initial  Condition  Constant,  K  for  Impact  Velocity,  4. A  km/s 


Parametric  Dependence  of  the  Initial  Condition  Conatant 


i 


Oi  (er* 


position 

e  M 


Figure  20.  Radial  and  Circumferential  Stress  Distributions  at  0.  1 
and  5  Microseconds  for  V  ■  3.3  km/s. 


after  the  penetrator  surface  passes  the  particular  axial  plane.  Both  radial 
and  circumferential  stresses  are  initially  compressive  near  the  cavity,  the 
circumferential  stress  is  less  than  the  radial  stress  by  the  amount, 

2 

(<r  +  M0/r  )  and  is  tensile  near  the  elastic/plastic  boundary.  As  time 
progresses,  both  stresses  decay  in  magnitude  and  at  5  microseconds,  the  radial 
stress  is  compressive  while  the  circumferential  stress  is  completely  tensile. 
The  spatial  gradients  of  both  stress  components  decay  very  quickly.  At  5 
microseconds,  the  stress  distributions  appear  almost  spatially  flat,  but  their 
difference  still  exceeds  the  static  yield  stress;  therefore  plastic  flow  will 
continue  until  the  yield  criterion  is  no  longer  met.  This  occurs  at 
approximately  23  microseconds  (Figure  19). 

The  selected  value  of  K  significantly  Influences  the  final  cavity 
dimension,  the  time  to  reach  this  value  and  the  Instantaneous  velocity  at  any 
radial  coordinate.  For  the  case  just  considered,  the  estimate  based  upon  the 
incompressible  flow  over  a  blunt  cylinder  yields  an  accurate  prediction  of  the 
final  cavity  dimension  for  both  Impact  velocities.  A  value  of  K  ~  .78  will, 
therefore,  be  assumed  for  the  remaining  parametric  calculations. 


B.  Dependence  upon  Impact  Velocity  and  Static  Target  Strength 


The  previous  section  presented  cavity  growth  solutions  to  equation  (92) 
for  the  case  of  a  tungsten  alloy  rod  penetrating  armor  plate.  The  solutions 
were  observed  to  depend  upon  the  initial  condition  parameter,  K,  which  relates 
the  initial  radial  velocity  with  the  steady  penetration  velocity,  U.  A  value 
of  K  ■  .78  produced  reasonable  agreement  between  the  calculation  and  the 
experiments  of  Silsby  [121].  Hohler  and  Stilp  [66]  performed  experiments  with 
both  steel  and  tungsten  alloy  rods  at  several  Impact  velocities  with  targets 
composed  of  different  strength  steels.  Figure  21  summarizes  their  results  for 
the  tungsten  rods  (1.5  ^  V  ^  4  km/s)  for  steels  with  hardnesses  of  180  BHN 

(or  ~  5  X  10^  dynes/cm^  and  300  BHN  (<r  ~  9  x  10^  dynes/cm^).  As  is 

y  y 

physically  expected,  the  normalized  final  cavity  dimension  decreases  with 
Increasing  target  hardness  (or  flow  stress>  and  increases  with  increasing 
Impact  velocity.  The  two  Isolated  points  represent  the  calculation  results 


from  the  previous  section  for  an  assumed  static  strength  of  (1  x  10 


10 


dynes/cm^)  and  a  slightly  larger  penetrator  dimension  (r  -  .39  cm).  The 


experimental  data  do  not  identify  the  dependence  of  penetrator  strength  and 
density  or  target  density,  but  the  calculations  asstjme  a  dependence  of  the 
form: 


y 


1 


1  + 


V 


so  as  the  penetrator  density  is  decreased,  the  penetration  rate  will  decrease 
and  one  should  then  observe  a  decrease  in  the  cavity  dimension  for  the  same 
Impact  velocity.  This  trend  is  observed  by  comparing  the  experimental  data 
presented  In  Figure  22  with  Figure  21.  The  results  for  the  high  hardness 
steel  (300  BHN)  show  negligible  penetrator  density  Influence  but  for  the  lower 
hardness  alloy  (180  BHN)  the  effect  is  more  pronounced  and  follows  the 
predicted  trend.  The  four  separate  points  on  Figure  22  represent  the  model's 
final  cavity  predictions  for  the  following  parameter  values: 
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Figure  21.  The  Normalized  Cavity  Dimension  Variation  with  Respect  to 
Impact  Velocity  and  Target  Strength  for  Tungsten  Rods  and 
Steel  Targets. 
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P  ■  7.8  g/cm"^  (measured) 

P  3 

P  -  7.8  g/cm  (measured) 

r  *  .27  cm  (measured) 

P  4 

/t  -  2  X  10  poise  (from  section  4-1) 

9  2 

<7  >  5  X  10  dynes/cm  (estimated) 

y  10  2 

to  1  X  10  dynes/cm  (estimated) 
V  «  3.5  X  10^  cm/s 

and  1.18  x  10^  cm/s 
K  -  .78. 


Once  again  the  initial  values  of  a  and  b  must  be  estimated  such  that  the 
Initial  radial  stress  at  a^  equals  the'^stagna^lon  pressure.  The  following 

values  were  used: 


V 

V 


3.5  km/s 

2.4  km/s 


9  2 

5  X  10  dynes/ cm 

8.8  X  10^  dynes/ cm^ 

9  2 

5  X  10  dynes/ cm 

1.0  X  10^^  dynes/ cm^ 


a 

o 

- 

.090 

cm 

b  -  .692 
0 

cm 

a 

.082 

cm 

b  -  .473 

cm 

o 

o 

a 

m 

.080 

cm 

b  -  .350 

rm 

o 

o 

a 

m 

.090 

cm 

b  -  .350 

cm 

o  o 


and  the  growth  curves  are  presented  In  Figure  23.  The  variation  of  the  final 
cavity  dimension  with  respect  to  both  the  static  yield  strength  and  Impact 
velocity  qualitatively  agrees  with  the  experimental  trend.  Quantitatively, 
the  agreement  Is  good  for  the  lower  strength  alloy  but  the  calculations 
overestimate  the  high  strength  alloy  In  Figure  22  and  underpredict  the  high 
strength  alloy  In  Figure  21.  Fart  of  this  discrepancy  Is  undoubtedly 
attributable  to  the  uncertainty  of  the  static  yield  strength.  Hohler  and 
Stllp  [66]  note  an  observed  range  of  measured  hardness  for  the  high  strength 
alloy  between  260  and  330  BHN.  The  quoted  yield  strength  varies  between  8.8  x 
9  10  2 

10  and  1.08  x  10  dynes/cm  .  The  growth  rate  model  Is  concerned  with  the 


flow  stress;  the  static  yield  strength  Is  only  one  component  of  this  property 
and  material  viscosity  Influences  the  Instantaneous  value.  Therefore,  the  two 
parameters,  and  m  must  be  more  precisely  defined  In  order  to  Improve  the 

agreement.  Elchelberger  [14]  presents  a  relationship  between  the  target 
Brlnell  hardness  and  an  "average  flow  stress"  defined  by  the  proportionality 
of  crater  volume  and  penetrator  kinetic  energy.  The  relationship  for  this 


g 

flow  stress,  v,  Is:  v  ■  2.6  x  10  x  (BHN).  Therefore,  the  expected  flow 
stress  for  the  high  strength  alloy  (BHN  >  300)  should  be  of  the  order,  8  x 
10  2 

10  dynes/cm  .  This  value  Is  an  order  of  magnitude  higher  than  the  static 


9  2 

yield  strength,  «r  »  9  x  10  dynes/cm  .  The  Hohenemser-Prager  yield 
assumption  related  the  flow  stress  with  the  static  yield  stress  In  the 


following  way: 


0  ■  Oy 


(147) 
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082  cm 


and  with  the  initial  values,  0(o}  1«12  x  10  cm/s  and  r(o)  »  a  >  . 

o 

for  the  high  strength  alloy  case.  As  shown  in  Figure  23,  the  influence  of  the 

9  6 

viscosity  coefficient  is  clear.  In  this  case,  <r  >  9  x  10  -)-;u(1.66x  10  ) 

2  4 

dynes/cm  .  A  value  of  m  *  2  x  10  was  determined  by  Harlow  [55]  and  used  in 
these  calculations.  Initially,  when  the  velocities  (rates  of  strain)  are 
greatest,  the  viscous  contributions  predominate  the  flow  stress.  The  sensi¬ 
tivity  of  the  calculations  to  the  viscosity  coefficient  will  be  Investigated 
in  the  next  section.  Leave  it  to  say  that  the  flow  stress  has  a  pronounced 
effect  upon  the  final  cavity  dimension  but,  unfortunately,  a  precise  value  for 
this  parameter  is  not  available  in  general. 

As  an  additional  check  of  the  overall  model,  calculations  of  cavity  growth 
in  aluminum  targets  are  presented.  Perez  [81]  published  the  results  of  Impact 
tests  involving  steel  rod  penetrators  into  2024  aluminum  targets  with 
velocities  ranging  between  2  and  4  km/s.  Perez  assumes  the  static  yield 

9  2 

strength  for  this  aluminum  alloy  as  2  x  10  dynes/ cm  .  The  Metals  Handbook 

9  2 

[124]  gives  the  tensile  strength  as  4  x  10  dynes/cm  ;  therefore,  both  values 
will  be  Independently  assumed  and  the  calculation  results  will  be  compared 
with  the  data  of  Perez.  Figure  24  presents  the  calculated  growth  curves  for 
the  following  problem  parameters: 

3 

P  «  7.8  g/cm 
P  3 

P  »  2.7  g/cm 

r_  ■  .15  cm 

P  q  2 

ff  •  2  X  10  dynes/cm 

^92 
4  X  10  dynes/cm 

V  ■  2.5  X  10^  cm/s 

4.0  X  10^  cm/s 

K  -  .78 

and  M  ■  5  X  10^  poise  (Walters  [125]). 


The  initial  condition  parameters  are: 

V  -  2.5  km/s  <T  ■  2  X  10^  dynes/cm^ 

^92 
(T  "  4  X  10  dynes/cm 

y  q  2 

V  ■  4.0  km/s  f,  ■  2  X  10  dynes/cm 

^92 
„  «  4  X  10  dynes/cm 


a 

m 

.050 

cm 

b 

m 

.281 

cm 

o 

o 

a 

o 

- 

.044 

cm 

b 

o 

m 

.186 

cm 

a 

o 

- 

.057 

cm 

b 

o 

- 

.551 

cm 

a 

0 

- 

.052 

cm 

b 

0 

- 

.340 

cm 

The  trends  are  similar  to  those  observed  in  the  steel  targets  with  the 
exception  that  the  times  required  for  cavity  arrest  are  consistently  shorter 
in  the  aluminum  target.  The  difference  is  most  likely  due  to  the  relative 
magnitudes  of  the  radial  momentum  distributions.  The  data  of  Perez  [81]  is 
reproduced  in  Figure  25  along  with  the  four  calculated  final  cavity 
dimensions.  The  model  predictions  again  qualitatively  agree  with  the 
experiment  but  since  the  flow  stress  (or  and  fi  parameters)  are  not 

precisely  known,  a  quantitative  comparison  should  not  be  attempted. 


As  noted  in  Appendix  B,  the  neglect  of  viscosity  effects  within  impact  and 
penetration  problems  is  not  Justifiable  in  general.  In  the  previous  section 
it  was  observed  that  the  flow  stress  has  an  important  role  in  determining  the 
final  cavity  dimension  and  that  both  the  static  yield  strength,  and  the 

target  viscosity,  Influence  the  specific  value  of  this  flow  stress.  This 
section  will  be  concerned  with  the  viscous  contribution  only.  No  experiments 
have  been  performed  that  Isolate  the  material  viscosity  parameter  and  its 
Influence  upon  the  target  cavity  dimension.  A  parametric  variation  of  the 
viscosity  coefficient  will  yield  a  variation  in  the  predicted  growth  behavior. 
A  particular  case  of  study  is  chosen  to  compare  with  the  finite  difference 
calculation  of  Kucher  [36]  and  the  problem  parameters  are: 

P  *8.9  g/cm^ 

P  3 

P  ■  7.8  g/cm 

r  -  0. 1  cm 

P  10  2 

ffy  ■  1.36  X  10  dynes/ cm 

V  ■  7.6  X  10^  cm/s 

K  -  .78 

0  2 

and  the  target  viscosities  considered  are,  in  poise,  m  >  1  x  10  ,  1  x  10  , 

1  X  10  and  1  x  10  . 

The  initial  conditions  for  the  first  three  viscosities  are  determined  from 
equation  (130)  and  for  the  last  viscosity  from  equation  (128)  since  the 
approximation  (129)  is  not  appropriate  for  this  order  of  magnitude.  The 
initial  values  are: 

M  -  1  X  10^,  1  X  10^,  1  X  10^  a  •  .036  cm  b  ■  .368  cm 

_  o  o 

»  ■  1  X  10^  a  »  .045  cm  b  ■  .461  cm 

o  o 

and  the  resulting  calculations  are  presented  along  with  the  calculation  of 
Kucher  [56]  in  Figure  26.^  The  effect  of  viscosity  is  not  significant  until 
values  of  the  order  of  10^  poise  or  greater  are  selected.  For  viscosities 

4 

greater  than  1  x  10  poise,  the  influence  is  not  linear.  An  order  of 
magnitude  variance  of  the  viscosity  coefficient  produces  a  reduction  of  only 
approximately  30Z  of  the  final  cavity  dimension.  Therefore,  reasonable 
estimates  should  be  expected  as  long  as  the  viscosity  coefficient  can  be 
determined  to  within  the  correct  order  of  magnitude.  Walters  [125]  notes  that 
given  the  current  state  of  understanding,  an  order  of  magnitude  estimate  is 
perhaps  as  good  as  one  might  expect.  The  relative  Importance  of  the  viscosity 
coefficient  can  be  observed  explicitly  in  the  relation  for  the  flow  stress, 
equation  (147).  Depending  on  the  magnitude  of  the  static  yield  strength,  the 
local  radial  velocity,  0/r  and  the  radial  position,  the  coefficient  may  or  may 
not  affect  this  flow  stress.  For  the  range  of  these  parameters  encountered  In 

this  example,  viscosity  on  the  order  of  10  poise  or  greater  will  contribute 
significantly  to  the  value  of  the  flow  stress. 


Figure  26.  Parametric  Dependence  of  the  Vlacoalty  Coefficient,  M  ,  Upon  the  Cavity  G.rowth  Behavior. 


The  finite  difference  calculation  of  Kucher  [56]  utilized  an  Eulerlan, 
compressible  and  explicitly  invlscld  formulation  but  the  computational  method 
introduced  an  Implicit,  efiEectlve  viscosity  of  undetermined  magnitude; 
although  It  Is  known  to  have  the  form 

A*  ■  l/2PUdx 

effective 

where  Ax  represents  the  characteristic  mesh  dimension  and  U  is  the  local  mean 
fluid  speed.  Assuming  the  results  of  Kucher  are  accurate,  the  value  of  m  ~  2 

4 

X  10  poise  (determined  by  Harlow  [55])  yields  the  growth  model  solution  with 
the  same  final  cavity  dimension.  However,  It  should  be  pointed  out  that  the 
two  differ  slightly  In  the  early  time  growth  rates. 


Dependence  upon  Penetrator  Radius 


The  governing  equations  (88,  92,  and  96)  Include  the  penetrator  dimension, 
r^,  as  a  problem  parameter.  One  would  expect  that  the  target  cavity  would 

Increase  If  the  penetrator  diameter  were  Increased  (all  other  parameters  held 
constant).  This  Is  reflected  In  the  normalized  data  published  by  Hohler  and 
Stllp  [66]  and  Perez  [81]  which  assert  a  linear  scaling  between  penetrator  and 
cavity  dimensions.  Frequently,  the  penetrator  dimension  Is  not  known 
precisely,  especially  for  shaped-charge  Jets.  For  these  situations  a  range  of 
values  Is  typically  assumed  and  predictions  are  similarly  determined.  As  an 
example  of  the  growth  model's  ability  to  handle  such  a  problem,  the  following 
parameters  are  assumed: 

Pp  -  8.9  g/cm^ 
p  ■  7.8  g/cm^ 

r  -  0.15,  0.2,  0.3  cm 
P  5 

V  ■  6.9  X  10^  cm/s 

K  -  0.78 

10  2 

V  «  1  X  10  dynes/cm 

^  4 

and  -  2  X  10  poise. 


This  case  represents  typical  high  speed  Jet  penetration.  Diperslo,  et  al. 
[126]  radiographically  measured  Jet  velocities  and  dimensions  at  selected 
times  after  charge  detonation.  The  Jet  was  copper  and  had  a  linear  velocity 
distribution  that  varied  between  8.3  and  1.0  km/s.  At  three  Inches  deep 
within  the  target,  the  Jet  was  eroded  to  an  Incoming  velocity  of  6.9  km/s  and 
the  Jet  radius  was  estimated  to  be  between  .15  and  0.3  cm  In  dimension.  The 
Initial  condition  parameters  can  be  Iteratively  determined  and  are: 


r  «  0.15  cm 
P 

r  -  0.2  cm 
P 

r  -  0.3  cm 
P 


a  >■  .056  cm 
o 

a  -  .074  cm 
o 

a  >  .112  cm 
o 


b 

o 

b 

o 

b 

o 


.601  cm 
.80  cm 
1.20  cm  . 


Figure  27  presents  the  calculated  growth  curves  for  the  specified  Input 
parameters.  The  measurements  [126]  were  taken  at  three  Inches  deep  Into  the 


Parametric  Dependence  of  the  Penetrator  Radius,  r  ,  upon  the  Cavity  Growth  Behavior 


target  (where  steady  conditions  would  be  likely)  and  represent  several  test 
shots  for  the  same  configuration.  The  cross  section  of  the  hole  at  each  axial 
position  is  not  perfectly  circular,  therefore  measurements  were  made  for  the 
greatest  and  the  least  diameter  at  each  axial  plane.  The  cavity  growth  model 
predictions  for  final  cavity  radius  identify  the  linear  dependence  upon  the 
penetrator  dimension.  The  time  required  to  reach  the  final  cavity  dimension 
also  exhibits  a  linear  dependence  upon  this  parameter.  The  experimental 
measurements  lie  within  the  range  of  the  calculated  predictions  for  the 
estimated  range  of  jet  radii. 

E.  Comparison  Between  Semi-infinite  and  Finite  Lateral  Boundaries 

The  preceding  sections  compare  postmortem  measurements  with  the  model 
predictions  for  ranges  of  the  problem  parameters.  All  cases  have  assumed  that 
the  target  dimension  in  the  plane  of  the  entrance  surface  is  large  in 
comparison  with  the  penetrator  diameter  or  cavity  dimension.  This 
configuration  is  typically  the  case  in  practice  but  occasionally  the  question 
arises;  if  the  lateral  dimension  is  not  overwhelmingly  large,  what  will  be  the 
effect  upon  the  growth  behavior  prediction?  The  formulation  of  the  governing 
equations  Includes  this  possibility.  Equation  (88)  was  derived  under  the 
presumption  that  the  lateral  boundary  is  circular  in  cross  section  and  is 
relatively  small  in  dimension.  The  radial  stress  is  assumed  to  vanish  at  this 
surface.  As  an  example  of  how  the  behavior  will  change  if  the  lateral 
boundary  is  small,  the  aluminum  target  case  of  section  IVB  will  be  recalled. 
The  following  parameters  are  assumed: 

P  ■  7.8  g/cm^ 

P  3 

p  ■  2.7  g/cm 

r  >  0.15  cm 
P  5 

V  *  4.0  X  10  cm/s 

K  -  0.78 

9  2 

9  ■  2  X  10  dynes/cm 

y  3 

and  /bi  -  5  X  10  poise. 


The  initial  condition  parameters  were  Iteratively  determined  for  the 
semi-infinite  case  from  equation  (134)  but  if  the  lateral  boundary  is  known 
and  the  radial  stress  is  set  to  zero  upon  this  surface,  then  equation  (128) 

must  be  utilized.  The  boundary,  b  ,  will  be  assumed  to  be  0.4  cm;  the  initial 

o 

cavity  position  is  then  iteratively  determined  to  be  0.050.  The  calculated 
result  for  this  finite  boundary  problem  is  presented  along  with  the 
semi-infinite  solution  in  Figure  28.  Initially  (t  <  5  ^s)  the  two  curves  are 
similar  but  as  time  proceeds,  a  difference  develops.  The  semi-infinite 
target  curve  exhibits  a  gradual  deceleration  of  the  cavity  until  the  motion  is 
arrested  due  to  the  diminishing  plastic  stresses.  Once  the  yield  criterion  is 
no  longer  met,  the  problem  is  terminated  since  only  elastic  displacements 
could  continue.  The  plastic  region  is  always  bounded  by  an  elastic  region  of 
infinite  mass. 

The  finite  lateral  boundary  problem  exhibits  noticably  different  behavior. 
The  deceleration  for  times  greater  than  5  microseconds  is  less  than  that  for 
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'v'-’w 


the  seml-lnflnlte  case.  It  is,  however,  not  zero  so  the  problem  will  continue 

for  a  much  longer  time.  At  40  microseconds,  the  difference  between  a  and  b 

o  o 

Is  .03  cm  at  a  radius  of  2.3  cm.  The  target  then  appears  as  a  thin,  expanding 
cylinder.  The  strain  experienced  by  the  tube  material  may  be  estimated  by 


T;  ^  -  46  at  t  -  40  MS. 

It  Is,  therefore,  likely  that  fracture  or  Instabilities  would  have  set  In 
by  this  time  In  an  actual  problem  and  that  the  original  plasticity  assumption 
may  not  even  apply.  It  would  be  difficult  to  perform  an  experimental 
measurement  of  the  final  "tube**  dimension  If  It  could  even  be  recovered 
Intact.  Perhaps  a  more  suitable  measurement  to  compare  with  this  calculation 
would  be  to  measure  the  outer  surface  velocity  In  an  average  manner  with  the 
use  of  a  high  speed  framing  camera.  At  this  point,  such  an  experiment  has  not 
been  performed.  From  the  calculation  results  alone.  It  can  be  qualitatively 
stated  that  decreasing  the  lateral  boundary  dimension  will  Increase  the  cavity 
dimension. 

F.  Comparison  with  Entrance  Hole  Data 

Throughout  this  chapter  the  calculated  final  cavity  dimensions  were 
compared  with  experimental  data.  Experiments  have  not  been  published  that 
measure  the  growth  behavior  during  the  expansion  process  deep  within  the 
target.  '  However,  radial  cavity  growth  at  the  Initial  Impact  surface  has  been 
observed  radiographically  [14,  35]  and  at  least  the  shapes  of  the  experimental 
curves  are  similar  to  the  shapes  of  the  predicted  curves  for  the  steady 
penetration  problem.  In  order  to  compare,  qualitatively,  the  Initial 
transient  growth  behavior  with  the  deep,  steady  growth  behavior,  the  following 
configuration  will  be  assumed.  Elchelberger  [14]  considered  the  entrance  hole 
expansion  In  an  aluminum  (HOOF)  target  block  resulting  from  the  Impact  of  a 
short  steel  disc,  V  >  5  km/s.  If  instead,  the  disc  thickness  Is  assumed  many 
times  longer,  then  a  steady  penetration  condition  should  occur  and  the 
equations  in  Section  IIF  should  then  apply.  The  experiment  Involving  the 
transient  entrance  surface  [14]  can  then  be  compared  with  the  predicted 
behavior  that  would  result  If  the  penetrator  had  greater  length.  The  problem 
parameters  are: 

3 

Pp  ■  7.8  g/cm 
P  -  2.7  g/cm^ 

r  *  0.4  cm 
P  5 

V  *■  5.0  X  10^  cm/s 

9  2 

>  2  X  10  dynes/cm  (assumed) 

K  >  .78  (assumed) 

3 

and  M  «  5  X  10  (as  In  Section  IVB). 

The  Initial  condition  parameters  are  Iteratively  determined  as: 

cm,  b^  ■>  1.928  cm  and  the  penetration  velocity,  U  Is  calculated  from  the 

steady  penetration  model  (Appendix  A)  as  3.15  x  10^  cm/s.  This  value  Is  less 
than  the  elastic  wave  speed  In  aluminum.  Therefore,  equation  (92),  the 
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subsonic,  seml-lnflnice  case,  will  be  utilized.  Figure  29  presents  both  the 
experimental  measurement  for  the  entrance  surface  and  the  calculated 
prediction  for  axial  planes  deep  within  the  target.  The  two  curves  are 
considerably  different  In  terms  of  both  final  dimension  and  Instantaneous 
slope  (velocity).  The  Isolated  sources  of  the  discrepancy  (Ignoring  the 
uncertainty  of  the  values  of  yield  strength  and  viscosity)  are  the  difference 
in  the  axial  boundary  conditions,  the  length  of  the  penetrator  and  the  degree 
of  transient  axial  behavior.  If  the  free  surface  boundary  condition  were  the 
only  variable  Involved,  It  would  be  expected  that  the  final  cavity  dimension 
for  the  steady  problem  be  less  than  the  final  dimension  for  the  entrance 
surface.  The  transient  penetration  rate  lies  between  the  Impact  velocity  and 
steady  penetration  rate;  consequently  the  Initial  condition  on  radial  velocity 
would  be  greater  and,  therefore.  It  would  be  expected  that  the  entrance  cavity 
be  larger  than  the  cavity  at  some  depth  Into  the  target.  However,  the  results 
show  Just  the  opposite  behavior.  This  Is  attributable  to  the  difference  In 
penetrator  length.  Elchelberger  fired  a  steel  disc  (radius  >  0.4  cm,  length  • 
0.05  cm)  end  on  at  the  seml-lnf Inlte  aluminum  target.  Steady  axial  behavior 
was  never  achieved  and  due  to  the  very  short  length,  relaxation  from  the  rear 
penetrator  surface  reduced  the  time  over  which  the  penetrator  could  Interact 
with  the  target.  Subsequently,  less  radial  Impulse  Is  Imparted  to  the 
entrance  plane  material  than  would  be  If  the  penetrator  length  were  Increased. 
Therefore,  In  order  to  study  the  effects  of  the  transient  nature  of  the  free 
surface  boundary  of  the  entrance  cavity,  the  length  of  the  penetrator  must  be 
Increased  and  the  experiment  must  be  repeated  (keeping  all  other  parameters 
fixed).  The  model  predictions  will  then  show  the  deviation  associated  with 
the  transient,  free  boundary  condition,  alone. 


V.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  RESEARCH 

The  primary  Intent  of  this  report  was  to  develop  a  model  for  the  expansion 
of  a  cavity  created  during  the  penetration  of  hypervelocity  rods  and  apply  It 
by  predicting  parametric  trends  for  problems  of  practical  Interest.  The  model 
was  formulated  by  making  several  assumptions  as  to  the  behavior  of  the  flowing 
target  material  during  the  cavity  expansion  process.  These  approximations  are 
supported  by  experimental  evidence  and  the  qualitative  results  of  finite 
difference  calculations.  By  specifying  stress  conditions  at  the  boundaries  of 
the  plastic  regions,  the  governing  system  of  nonlinear  field  equations  was 
reduced  to  simpler  ordinary  differential  equations.  Upon  estimation  of  the 
Initial  conditions,  the  Initial  value  problems  were  solved  numerically  for 
several  situations  where  experimental  results  were  available.  The  following 
conclusions  can  therefore  be  stated: 

a.  The  calculated  predictions  of  the  final  radial  cavity  dimension  show 
the  same  parametric  variation  as  the  experiments. 

b.  The  time  dependent  growth  of  the  cavity  agrees  qualitatively  well  with 
the  finite  difference  calculations  of  Kucher  [56]  but  Is  significantly 
different  from  the  experimentally  observed  entrance  hole  behavior  of  short 
disc  Impact  [14]. 


c.  The  accuracy  of  the  model  predictions  depends  upon  the  accuracy  with 
which  the  initial  conditions  (V,  r^,  K,  U)  and  the  target  material 

characteristics  fi)  may  be  determined. 

d.  The  analogy  with  the  incompressible  fluid  flow  over  a  blunt  cylinder 
yields  an  estimation  for  the  initial  conditions  that  consistently  results  in 
close  agreement  with  the  experiments. 

e.  The  Hohenemser-Prager  vlsco-plastlc  f low  assumption  predicts  initial 
flow  stress  values  within  the  range  of  those  published  by  Elchelberger  [14] 
from  Independent  experiments. 

f.  The  analysis  is  applicable  to  cavity  growth  resulting  from  both 
kinetic  energy  penetrators  (1  km/s  ^  V  ^  4  km/s)  and  shaped  charge  jets  (2 
km/s  <.  V  ^  10  km/s)  as  long  as  the  physical  characteristics  of  each  is  known 
to  sufficient  accuracy  and  the  steady  penetration  presumption  is  valid. 

The  formulation  required  axial  symmetry,  material  Isotropy,  plastic 
incompressibility  and  steady  penetration.  If  any  of  these  approximations  is 
violated,  a  new  formulation  must  be  performed.  As  already  noted  in  Appendix 
A,  the  steady  penetration  theory  may  not  be  applicable  to  material 
combinations  with  substantially  different  compressibilities  or  when  the 
hydrodynamic  pressure  is  of  the  order  of  the  target  strength  or  less. 
Similarly,  the  radial  cavity  expansion  model  is  limited  to  problems  where 

9 

hydrodynamic  penetration  is  applicable  (V  >  1  km/s  for  a  ~  1  x  10 

2  ^ 
dynes/ cm  ).  The  model  is  limited  also  by  the  requirement  for  input  values 

that  may  not  be  easily  determined,  i.e.,  m,  (Xy,  etc.  Additionally, 

materials  that  do  not  exhibit  plastic  flow  prior  to  rupture  or  fracture  cannot 
be  investigated  since  the  plastic  characteristics  have  implicitly  been  assumed 
to  predominate  the  cavity  expansion  process. 


.  Beyond  the  obvious  requirement  for  further  research  of  material 
characterization  during  ballistic  Impact,  an  alternate  formulation  of  the 
completely  transient,  initial  Impact  problem  for  long  rod  penetrators  should 
be  pursued  and  be  compared  with  the  results  from  this  report.  An  experimental 
observation  of  the  growth  of  the  cavity,  deep  within  the  target,  or  stress 
measurements  (similar  to  those  of  Pritchard  [70]}  at  points  displaced  from  the 
penetration  trajectory,  taken  during  the  early  expansion  times,  are  also  very 
much  needed. 
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LIST  OF  SYMBOLS 

(Units  According  to  COS  Convention) 
radial  position  of  the  cavity  surface 
initial  value  of  a 

radial  position  of  the  boundary  that  separates  Regions  II  and  III 
radius  of  penetrator  bulb,  equation  (lOO) 

Bethe's  flow  stress  parameter,  equation  (105) 
transform  variable  '.n  Section  (IIIB) 
cross  sectional  areas 

radial  position  of  the  elastic/plastic  boundary 

initial  value  of  b 

body  force  vector 

Bethe's  radial  pressure  parameter 

radial  position  of  the  target  lateral  dimension 

wave  speed 

elastic  wave  speed 

integration  constant  in  Section  (IIIB) 

Incremental  mechanical  work 
incremental  cavity  volume 

Incremental  control  volume  deflection  angle 

body  of  revolution  diameter 

cavity  diameter 

penetrator  diameter 

strain  rate  tensor 

elastic  component  of  strain  rate  tensor 
plastic  component  of  strain  rate  tensor 
yield  function 


target  resistance  force 

initial  condition  parameter,  equation  (ll6) 

elastic  boundary  constant 

elastic  boundary  constant 

transform  constant  in  Section  (IIIB) 

stagnation  velocity  gradient 

maximijm  shear  rate 

length  of  penetrator  bulb 

transform  variable  in  Section  (IIIB) 

penetrator  mass 

local  Mach  number 

transform  variable  in  Section  (IIIB) 

pressure  in  the  penetrator  control  volume 

transform  variable  in  Section  (IIIB) 

stagnation  pres  sure 

penetration  depth  in  Appendix  A 

radial  coordinate 

penetrator  radius 

radius  of  curvature 

relative  velocity  scalar 

also  root  to  principle  stress  determinant 

relative  velocity  vector 

principle  stress  vector 

•  penetration  rate 

velocity  at  the  edge  of  the  boundary  layer 

•  velocity  vector 


particle  velocity 


local  mean  fluid  speed 
displacement  vector 
impact  velocity 

distance  from  stagnation  point 
coordinate  vector 

transform  variable  in  Section  (IIIB) 
axial  coordinate 

transform  variable  in  Section  (IIIB) 

shear  strain  in  the  r-z  plane 

Kronecker  delta 

characteristic  mesh  dimension 

circumferential  coordinate 

angular  displacement  from  symmetry  axis 

deflection  angle  of  microstructure  bands 

slope  of  the  U  -  data 
3  p 

bulk  viscosity  coefficient,  equation  (U6) 

empirical  fitting  parameter,  equation  (5**) 

shear  viscosity  coefficient 

Poisson's  ratio 

3.1U  .  .  . 

target  density 

penetrator  density 

target  density 

average  flow  stress  during  high  strain,  strain  rate  deformation 

target  yield  strength  in  uniaxial  tension,  measured 
quasi-statically 


stress  tensor 
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APPENDIX  A 


DERIVATION  OF  THE  STEADY  PENETRATION  RATE 

The  theory  describing  the  penetration  of  high  speed  rods  and  Jets  was 
developed  during  World  War  II  [58,  79.  127  -  129]  presumably  to  allow  the 
design  of  armor  and  anti-armor  warheads*  Working  independently.  Hill,  Mott 
and  Pack  [79]  (England.  1944)  and  Pu^  [127]  (U.S.A.,  1944)  considered 
essentially  the  same  configuration,  as  shown  in  Figure  A7, 


Figure  Al.  The  Penetration  of  a  Semi-Infinite  Target  with  a  Long  Penetrator. 


where  Ap  represents  the  penetrator  cross  sectional  p  and  p  p  the  taraet  and 
penetrator  densities,  respectively.  V  represents  the  velocity  of  the 
penetrator  with  respect  to  a  stationary  frame,  A  represents  the  crossectlon  of 
the  target  over  which  the  penetrator  interacts  and  the  surface  W  -  W' 
represents  a  shock  surface  (if  the  velocity  is  sufficient).  The  interface  is 
presumed  to  translate  to  the  right  at  a  constant  velocity,  (J.  If  the  frame  of 
reference  is  attached  to  the  Intersection  point  of  penetrator  surface  and  the 
symmetry  axis,  moving  with  a  constant  speed  U,  then  the  axial  flow  in  this 
moving  frame  appears  steady.  Surface  A  is  bombarded  on  one  side  by  the 
penetrator  moving  with  a  relative  velocity  (V  -  U)  and  on  the  other  aide  by 
target  material  moving  in  the  opposite  direction  with  relative  velocity  U. 
Within  this  system  of  translating  coordinates,  the  incoming  penetrator  is 
sprayed  out  radially  upon  impact,  losing  its  axial  momentum.  Therefore,  the 
forces  on  both  sides  of  A  must  be  equal  and  opposite.  Pugh,  thereby  obtained: 


PpAp(V  -  U)^  =  PAU^ 


(A.1) 


and  upon  solving  for  U, 


Equation  (A. 2)  can  be  further  alnpllfied  by  aeauming  that  A  s  kp,  then 


(A. 3) 


Hill  obtained  the  same  relationship  for  the  penetration  velocity  by 
applying  the  theory  of  free  streamlines.  Upon  the  translating  stagnation 
point,  the  pressures  (given  by  the  Bernoulli  equation)  must  equal  In  both  the 
target  and  penetrator.  For  the  lneoiiq)resslble  ease. 


1/2Pp(V-0)^  s  1/2pU^ 


(A. 4) 


and  equation  A. 3  then  follows.  Both  Hill  and  Pugh  recognized  that 
compressibility  of  either  the  penetrator  or  target  would  Influence  the  value 
of  the  pressure  on  either  side  of  the  stagnation  point,  i.e.. 


(A. 5) 


The  effects  of  compressibility  were  subsequently  studied  by  Allison  and 
Vltall  [130],  Harlow  and  Pracht  [551,  Haugstad  and  Dullum  [131,  132]  and  Fils 
and  Chou  [133],  Considering  metal  targets,  Allison  concludes  that  the 
Incompressible  approximation  Is  sufficiently  accurate  for  practical  problems. 
Harlow,  using  a  stlffened-gas  equation  of  state  for  steel  and  aluminum  with  a 
multidimensional  Eulerlan,  finite  difference  technique,  concludes  that  both 
the  steady  and  Incompressible  assumptions  are  quite  appropriate.  The  works  of 


Haugstad  and  Flis,  however,  ahow  that  for  materials  with  significantly 
different  compressibilities,  the  error  resulting  from  the  application  of  the 
incompressible  assun^tlon  to  the  determination  of  the  pressure  and  hence 
penetration  velocity  could  be  on  the  order  of  30Z. 

The  presuiiq)tion  that  the  materials  behave  like  fluids  is  based  on  the 
argument  that  the  inherent  strength  of  the  target  and  penetrator  materials  is 
much  less  than  the  hydrodynamic  pressure.  For  lesser  Impact  velocities,  this 
presumption  is  not  valid  and  equation  (A.t)  is  usually  modified  to  include 
strength-like  terms  and  a  pressure  variation  term, y  , 


pVi^  aYp  (V  -  0)^  +  * 


(A. 6) 


With  an  estimate  of  the  "resistance  stress,"  <r  and  a  guess  as  to  the  value 
ofy,  1<y^2  (Eichelberger  [62],  Allen  and  Rogers  [134])  the  penetration 
velocity  can  be  determined  once  again.  The  depth  at  which  the  penetration 
ceases  can  be  estimated  by  assuming  that  the  process  ends  idien  the  trailing 
end  of  the  penetrator  reaches  the  stagnation  point.  Integrating  the 
penetration  velocity  with  respect  to  time  yields  the  depth. 


r 

P(t)p  *  j  0  dt»  s  ^  yp  I 


(V  -  0)dt  s 


(A. 7) 


where  L  represents  the  constant  penetrator  length.  Eichelberger  [62] 
performed  an  experimental  test  of  this  theory  for  the  case  of  metallic  Jets. 
The  major  distinction  between  a  Jet  and  a  hypervelocity  rod  is  that  the  Jet 
has  a  variable  velocity  along  its  length.  This  casts  doubt  on  the  steady 
approximation  and  it  would  appear  that  the  above  equations  might  not  be 
applicable  to  these  types  of  penetrators.  However,  Eichelberger  concludes 
that  the  hydrodynamic  approximations  are  sui*prlslngly  accurate  when  compared 
with  experimental  evidence.  This  accuracy  should  diminish  with  decreasing 
impact  velocity  where  the  solid  strength  characteristics  predominate. 

However,  Christman  and  Gehring  [24]  compared  the  steady  penetration  theory 
with  experimental  data  for  high  speed  (1-6  km/s)  rods  and  concluded  that 
within  the  range  of  variables  considered,  the  target  strength  had  negligible 
influence  on  the  penetration  rate.  This  is  also  confirmed  by  the  finite 
difference  calculation  of  Van  Thiel  [135].  In  review  of  the  many  experimen¬ 
tal  observations  and  analytic  predictions,  it  can  be  assumed  that  the  steady, 
incompressible  penetration  model  produces  a  sufficiently  accurate  estimate  of 
the  penetration  rate  of  long  penetrators  into  semi-infinite  metallic  targets. 


APPENDIX  B 


THE  YIELD  STRESS  DURING  BALLISTIC  IMPACT 


Perhaps  one  of  Che  firac  researchers  to  be  concerned  vic^'  the  effect  of 
the  yield  stress  variation  upon  the  ability  of  a  penetrator  to  perforate  armor 
was  H.  A.  BeChe  [19]  who  published  in  1941.  In  the  formulation  of  his  armor 
penetration  model,  he  introduced  the  Tresca/Mohr  yield  criterion,  which 
supposes  that  the  target  material  exhibits  an  observable  elastic-plastic 
transition  stress,  cr^.  This  yield  stress -determined  the  magnitudes  of  the 

resistance  stresses  within  the  armor  and  consequently  became  an  important 
problem  parameter.  Bethe  recognized  that  previous  theoretical  works,  (Becker 
[136]  who  published  in  1926  and  Orowan  [137]  who  published  in  1935)  although 
not  nearly  complete,  had  predicted  through  statistical  mechanics  arguments 
chat  Che  yield  stress  at  high  rates  of  strain  should  be  much  greater  Chan  Che 
yield  stress  observed  in  quasi-static  testing.  In  this  pioneering  paper, 

Bethe  included  a  review  of  the  early  experimental  observations  of  Hopkinson 
[138],  Elam  [139],  Davis  [140]  and  Manjoine  and  Nadai  [141]  which  all 
confirmed,  at  least  qualitatively,  the  aforementioned  theoretical  trend.  The 
experiments  of  Elam  and  Davis  showed  that  even  chough  Che  yield  stress 
increased  up  to  about  one  hundred  sixty  percent  of  the  static  value  (limited 
by  Che  speed  of  the  testing  machine),  it  did  so  with  practically  no  strain 
hardening.  The  observations  allowed  Bethe  to  select  an  analytically  tractable 
form  for  subsequent  calculations.  Interest  in  discerning  the  behavior  of 
materials  undergoing  sianiltaneous  large  strains  and  rates  of  strain  seemed  to 
escalate  during  the  years  following;  many  experiawntal  studies  were  published 
and  many  constitutive  equations  were  proposed.  Reviews  of  these  works  can  be 
found  in  several  current  texts  (CrisCescu  [88],  Billington  and  Tate  [2], 
Johnson  [4],  Zukas,  eC.al.  [6]).  The  conclusion  reached  in  all  instances  is 
chat  the  dynamic  yield  stress  is  not  a  characteristic  constant  of  the 
material,  its  value  depends  upon  the  rate  of  loading  as  well  as  the  form  of 
Che  constitutive  equation  used  for  a  considered  material.  Quoting  Cristescu, 
"the  mechanism  by  which  the  loading  rate  influences  the  static  stress-strBin 
curve  is  not  yet  clearly  understood."  This  is  still  true  today  and  can  be 
observed  by  comparing  the  research  papers  of  the  rate-independent,  strain 
hardening  proponents  (Bell  [142,  143])  with  chose  of  the  rate-dependent 
community  (Malvern  [144],  Perzyna  [145],  Gillis  [146],  Klahn  [147],  etc.  and 
the  numerous  Russian  workers  reviewed  by  Walters  [125]).  From  dynamic  wave 
propagation  experiments.  Bell  proposes  that  the  following  stress-strain, 
function  has  "wide  and  remarkable  generality"  for  strain  rates  up  to  10^ 

---“I 


(”r) 


G(o)  B 


(B.l) 


where  G(0)  is  the  static,  isotropic,  linear  elastic,  shear  modulus,  r*  is  an 

integer  index,  and  B^  is  a  universal  dimensionless  constant.  The  rate  of 

o 

strain  does  not  appear  explicitly  within  this  expression.  This  is  in 
contradiction  to  the  many  dynamic  yield  .experiments  performed  by  many 
researchers  over  many  years  (Alder  and  Phillips  [148],  Butcher  and  Karnes 
[149],  Campbell  [150,  151],  Hawkyard  [152],  Hecker  [153],  Malvern  [144], 
Taylor  [154],  Whiffen  [155],  Wilkins  [156],  White  [157],  Zener  [158]  and  many 
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others).  Alder  and  Phillips  [148]  perforned  experiments  oh  aluminum,  steel 

and  copper  with  strain  rates  up  to  40  sec~  and  found  they  could  best  fit  the 
data  with  a  power  law  form: 


(B.2) 


where and  n  are  experimentally  determined  constants.  Ludwlk  [159]  proposed 


the  following  logarithmic  equation: 


o  in 


(B.3) 


where  <r^  and  are  constants,  and  the  plastic  strain  rate  Is  normalized  with 

respect  to  the  quasl-statlc  strain  rate.  Malvern  [144]  extended  this  to 
Include  work-hardening,  <r  s  f(e). 


9  *  fCe)  ♦  ♦  be^). 
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Hohenemser  and  Prager  [107],  and  later  Kotllarevskll  [160]  and* Perzyna  [145] 
assumed  that  rate  sensitive  plastic  materials  could  be  modeled  with  a  Bingham 
type  constitutive  form. 


(B.5) 


where  the  first  term  Is  the  rate  Independent  yield  stress  and  that  the  second 
term  represents  the  viscous  contribution  of  the  yield  stress  dependent  on  the 
material  viscosity,  which  most  likely  varies  with  the  current  state  of  the 
material.  Many  Russian  researchers  (see  Walters  [125])  have  utilizled  this 
last  equation  for  the  yield  stress  dependence  on  strain  rate.  Godunov,  et. 
al.  [105]  reviews  both  the  experimental  evidence  and  the  Russian  modeling 
philosophy  for  many  materials  Including  the  typical  metals.  Godunov  [106],  In 
a  later  report,  notes  the  variation  of  the  experimentally  determined  viscosity 
coefficient  with  different  rates  of  deformation  (strain  rates)  and  concludes, 
as  In  Figure  Bl,  that  the  viscosity  coefficient  has  only  weak  dependence  on 

.  ii 

the  rate  of  strain  If  the  strain  rate  Is  sufficiently  great,  l.e.  e  >  10  . 

This  then  simplifies  equation  (B.5)  considerably.  If  the  viscosity 
coefficient  can  be  Independently  determined,  given  of  course  the  expected 
temperature  and  pressure  regime,  then  the  dynamic  yield  stress  (or  flow  stress 
as  It  Is  alternately  referred  to)  is  a  linear  function  of  the  rate  of  strain 
for  a  specific  range  of  this  rate.  The  coefficient  of  viscosity,  being 
approximately  constant,  would  change  when  passing  from  one  range  to  another. 
Perzyna  [145],  recognizing  the  Insignificance  of  strain  hardening  In  mild 
steel  (Campbell  [150,  151]),  generalized  the  constitutive  equations  of 
Hohenemser  and  Prager  [107]  to  general  states  of  stress. 


For  F  >  0 


‘ii* 

6  - 

*ij  -  IT- 


For  F  <  0 


where  F  is  the 
deviation,  e 

•  J 

and  R,  as  well 


(B.6) 
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static  yield  function,  denote  the  components  of  stress 
the  components  of  strain  deviation  and  the  elastic  constants,  G 
as  the  function  y4>(F)  must  be  determined  experimentally. 


The  Influence  of  viscosity  upon  the  impact  and  penetration  problem  has 
been  Investigated  theoretically  by  several  authors.  Rlney  [51,  52], 

MacCormack  [161],  Walsh  [162]  using  multidimensional,  explicit,  finite 
difference  techniques  concluded  that  viscous  effects  are  quite  significant 
especially  when  characteristic  dimensions  are  small.  The  analytic  models  of 
Selzew  [43],  Whitesides  [47]  and  Ravld  and  Bodner  [163]  recognize  that  viscous 
effects  predominate  the  flow  In  the  later  stages  of  cavity  development. 

During  the  later  stages,  the  cavity  approaches  Its  final  state  and  the  effects 
of  the  Initial  hydrodynamic  Impact  stage  are  minimal.  The  neglect  of  rate 
dependence  Is,  therefore,  difficult  to  justify  In  general.  The  National 
Materials  Advisory  Board  of  the  National  Academy  of  Sciences  published  a 
summary  report  [164]  that  resulted  from  an  extensive  review  of  current  consti¬ 
tutive  modeling  efforts  oriented  towards  utilization  within  the  modern 
numerical  finite  difference  methods.  It  was  concluded  that  the  sophisticated 
numerical  methods  are  ultimately  limited  by  the  current  uncertainty  In 
material  response  descriptions.  A  relatively  poor  understanding  of  the 
mlcromechanlcal  mechanisms  of  dynamic  plastic  deformation  forces  the  use  of 
simple,  phenomena logical  continuum  models.  The  assumption  of  perfect 
plasticity  with  a  constant  flow  stress  has  been  found  to  yield  "excellent 
agreement  for  a  number  of  ordnance  designs  and  is  almost  always  used  In 
ordnance  calculations  at  present,"  even  though  the  actual  problems  are 
characterized  by  geometric  and  material  complexities.  The  flow  stress  Is 
always  different  from  the  statically  determined  value,  but  Its  precise  value 
Is  not  well  established.  The  flow  stress  may  have  explicit  dependence  on 
temperature,  strain,  strain  rate,  or  pressure  in  a  variety  of  functional 
forms. 


Equations  (B.1  -  B.6)  represent  only  a  partial  sampling  of  the 
constitutive  equations  that  have  been  assumed  in  order  to  solve  actual 
problems  Involving  Impact  and  the  resulting  plastic  deformation.  The 
behaviors  of  specific  materials  are  often  best  modeled  by  only  one  of  the 
above  equations,  and  even  then,  only  after  Judicious  choice  of  the 
undetermined  constants.  The  ultimate  decision  of  which  constitutive  equation 
to  use  Is  usually  based  on  considerations  of  analytic  tractablllty, 
reproduction  of  experimental  evidence  and  compatibility  with  theoretical 
arguments.  The  material  constants  must  also  be  available  If  quantitative 
comparisons  are  desired  In  addition  to  the  parametric  dependencies.  It  should 
be  noted  at  this  point  that  errors  resulting  from  the  approximation  of  the 
material  behavior  will  moat  likely  manifest  themselves  as  dlscrepencles 
between  analytic  predictions  and  experimental  observations.  With  this  in 
mind,  review  of  past  modeling  efforts  and  their  respective  successes  should  be 
completed  before  selection  of  a  particular  constitutive  model  for  a  specific 
problem. 
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APPENDIX  C 


THE  SOLUTION  OF  THE  AXIAL  VELOCITY  DISTRIBUTION 

The  formulation  of  Region  II  Involved  approximations  to  the  three 
dimensional  flow  field  that  were  both  experimentally  and  numerically 
supported.  Equation  (65)  originates  from  the  axial  momentum  equation  and 
couples  the  temporal  and  spatial  behavior  of  the  axial  velocity  component  as: 


3“z  .la^j  M 
at  •  7  ar  j  p 


(C.1) 


where  0(t)  Is  the  ultimate  solution  function  determined  by  solving  the  Initial 

value  problems  of  Section  II. F.  Suppose  that  the  axial  velocity  can  be 
represented  by: 


U^  =  T(t)R(r) 


(C.2) 


where  T  and  R  represent  functions  of  time  and  radial  coordinate  respectively. 
Then  equation  (C.1)  becomes: 


51  -  1  SSi  ii.  . 
at  *  r  ar  I  7  " 


(C.3) 


Following  arguments  standard  to  the  method  of  separation  of  variables, 


T(^-  0(t)) 


21 

dt 


-2—  —  -  constant  s  -  <1/ 
Rr  dr 


(C.4) 


the  partial  differential  equation  (C.1)  has  been  reduced  to  two  ordinary 
differential  equations  (C.4).  Integrating . both  of  these  equations  (C.4) 
yields 


and 


0(t) 


R  s  exp 


♦  V 


(C.5) 


(C.6) 
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Therefore,  via  equation  (C.2)  the  axial  velocity  may  be  determined: 


0 

z 


A 

-  f  {JLt  -  0{t) 

~2 


(C.7) 


0(t)  must  be  independently  determined,  is  determined  by  the  application 
of  the  initial  condition. 


U  (r  =  a  ,  0)  =  0 
z  o 

0(0)  =  0 

and  equation  (C.7)  then  becomes: 


*  0  exp  j  -  0(t)  ♦  r^^  ^  (C.8) 


The  eigen  value  f  remains  undetermined  but  could  be  empirically  deduced  from 
the  experimental  data  presented  in  Section  II. D. 
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APPENDIX  D 


DESCRIPTION  OF  THE  SOLUTION  PROCEDURE 


Determine  the  problem  parameters, 

‘'p'  '  ^p'  ^o'  ^ 

Calculate  the  penetration  rate,  U,  from  equation  (A. 3) 

Determine  the  appropriate  case:  supersonic  penetration,  subsonic 
penetration  and  aeml-lnflnlte  lateral  dimension  or  subsonic  penetration 
and  finite  lateral  dimensions. 

Select  a  value  of  K  (K  =  ^4  Is  suggested)  and  determine  the  remaining 
Initial  values,  a^  and  b^  by  the  procedure  outlined  Infection  II. G  (for 

example:  equation  (134)). 

5.  Select  the  appropriate  governing  equation  presented  In  Section  II. F: 
supersonic,  (96) 


1+  U^(a^+20)\ 

.-f-) 


subsonic  -  finite  lateral  dimension 


(b:+22) 


(b%20) 


0(a  ♦20) 


0-’f1+  U‘(a  >20) 


subsonic  -  seml-inflnlte  lateral  dimension,  (92) 


The  equation  (92)  la  Included  In  the  FORTRAN  program  In  Table  1.  The 
other  two  cases  can  be  considered  with  the  programming  modifications 
outlined  in  Section  ill, a.  * 

6.  Input  all  the  necessary  terms.  In  CCS  units,  and  choose  an  Initial  time 

Increment  of  10*^^  seconds.  After  the  calculation  Is  complete,  check  the 
solution's  first  cycle  such  that  the  Initial* pressure  Is  that  of  the 
stagnation  approximation.  Adjust  the  time  increment  to  confirm 
convergence. 
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